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ABSTRACT 


This  thesis  takes  five  FORTRAN  IV  programs  from  MCom- 


puter  Programs  for  Computational  Assistance  in  the  Study 


of  Linear  Control  Theory’’  by  Melsa  and  Jonc^p  and  translates 
them  into  a  microcomputer  BASIC  language  to  run  on  an  inex¬ 
pensive  microcomputer  system.  Three  of  the  five  programs  are 
state  variable  programs.  They  are  BASMAT  for  basic  matrix 
manipulation,  RTRESP  for  rational  time  response,  and  GTRESP 
for  graphical  time  response.  Two  are  transfer  function  pro¬ 
grams,  FRESP  for  frequency  response  and  RTLOC  for  root  locus. 
A  user's  guide  and  example  are  included  for  each. 

A  final  example  is  used  to  demonstrate  the  utility  of 


the  two  transfer  function  programs  as  an  aid  to  direct 
digital  design  in  the  w' -plane. 
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I .  INTRODUCTION 


The  purpose  of  this  thesis  is  two  fold.  First,  it 
is  to  show  that  standard  computer  programs  useful  in  the  study 
of  linear  control  theory  may  be  adapted  to  run  on  an  inex¬ 
pensive  home /microcomputer .  To  demonstrate  this,  five 
programs  were  chosen  from  [Ref.  1],  three  state  variable 
programs  and  two  transfer  function  programs.  The  three  state 
variable  programs  are  BASMAT,  a  basic  matrix  manipulation 
program,  RTRESP,  a  rational  time  response  program,  and  GTRESP, 
a  graphical  time  response  program.  The  two  transfer  function 
programs  are  FRESP,  a  frequency  response  program,  and  RTLOC , 
a  root  locus  program.  These  programs  as  they  appear  in  [Ref. 
1]  are  written  in  basic  FORTRAN  IV  language  to  be  run  on 
a  main  frame  computer  system  utilizing  standard  graphics 
subroutines.  These  programs  were  modified  and  rewritten  in 
a  microcomputer  BASIC  language  which  is  an  interpreted 
language.  Generally  these  programs  are  limited  to  systems 
of  10th  order.  It  is  felt  that  this  limitation  is  more  than 
acceptable  for  the  purposes  of  this  thesis.  Also  no  major 
effort  has  been  made  to  analyze  or  improve  the  efficiency  of 
the  numerical  methods  used.  An  effort  of  this  type  is 
advised  if  these  programs  are  to  be  modified  for  higher 
order  systems. 


The  second  part  of  this  thesis  investigates  the  methods 
and  relationships  involved  in  direct  digital  design  in  the 
w'-plane.  The  two  transfer  function  programs  adapted  to  run 
on  a  microcomputer  from  the  first  part  of  this  thesis,  FRESP 
and  RTLOC,  are  used  to  aid  in  this  investigation. 

Section  II  summarizes  the  common  problems  and  consi¬ 
derations  involved  in  the  translation  of  the  five  programs 
to  be  run  on  a  microcomputer  system.  Section  III  describes 
the  three  state  variable  programs,  BASMAT,  RTRESP,  and 
GTRESP,  and  gives  an  example  of  their  use  and  output. 

Section  IV  is  similar  to  section  III  and  describes  the  two 
transfer  function  programs,  FRESP  and  RTLOC. 

Section  V  deals  with  the  w'-plane.  Subsection  A  gives 
some  background  on  the  w'-plane  and  subsection  B  provides  a 
simple  example  using  the  two  transfer  function  programs  to 
compare  the  s  and  w'  planes.  Subsection  C  develops  templates 
of  some  constant  parameters  in  the  w'-plane.  Subsection  D 
ends  the  section  with  a  more  involved  example. 

Section  VI  provides  some  conclusions  and  recommendations. 


TRANSLATION  CONS IDERAT I ONS 


The  five  programs  in  this  thesis,  BASMAT,  RTRESP,  GTRESP, 
FRESP,  and  RTLOC ,  were  translated  directly  from  programs  of 
the  same  name  in  [Ref.  1].  These  programs  were  written  in 
the  basic  FORTRAN  IV  computer  language,  a  compiled  language 
designed  to  run  on  a  main  frame  or  mini  computer  system.  In 
the  first  part  of  this  thesis  these  programs  are  translated 
into  a  microcomputer  BASIC  language  adapted  to  run  on  the 
microcomputer  system  described  in  Appendix  A. 

In  general  the  translation  from  FORTRAN  IV  to  the  micro¬ 
computer  BASIC  used  posed  no  serious  problems.  Even  though 
the  BASIC  language  used  is  by  necessity  an  abridged  version 
of  the  BASIC  language  found  on  larger  more  expensive  systems, 
it  was  extensive  enough  to  provide  the  necessary  commands 
for  these  programs. 

The  input  portion  of  the  programs  were  rewritten  to  be 
interactive  for  convenient  keyboard  entry  of  problem  para¬ 
meters  eliminating  input  formatting  errors.  The  programs 
requiring  more  extensive  input  were  given  the  added  capability 
of  saving  and  retrieving  a  problem  description  from  a  disk 


file  for  the  user's  convenience. 


Unlike  FORTRAN  IV  the  BASIC  used  does  not  have  complex 
math  capabilities.  Therefore,  mathematical  operations  on 


complex  quantities  were  programed  separately  for  the  real  and 
imaginary  parts. 


Another  important  difference  in  the  languages  affecting 
translation  is  that  BASIC  has  no  provision  for  local  variables. 

All  variables  in  the  main  program  and  all  subroutines  are 
global.  This  caused  some  bookkeeping  problems  in  translating 
subroutines  to  prevent  undesired  side  effects.  Another  compli¬ 
cation  was  that  FORTRAN  IV  considers  the  first  four  characters 
of  a  variable  name  for  identification  of  that  unique  variable. 

The  BASIC  used  considers  only  the  first  two.  This  further 
complicated  the  bookkeeping  of  variables  and  resulted  in 
variable  names  being  assigned  just  because  they  were  different 
from  the  rest  and  with  no  relation  to  the  quantity  represented. 

The  output  routines  were  written  to  conform  as  close  as 
possible  to  the  FORTRAN  IV  version.  Since  the  microcomputer 
system  used  is  limited  to  eighty  columns  on  eight  and  a  half 
inch  wide  paper,  provisions  were  made  within  the  program  to 

4 

automatically  switch  to  a  condensed  character  font  when 
necessary  to  output  greater  than  eighty  characters  per  line. 

The  programs  requiring  graphical  output  created  some 
unique  problems.  First,  no  standard  library  subroutines  for 
graphics  as  used  in  the  FORTRAN  IV  programs  were  available 
on  the  microcomputer  system  used.  As  a  result  all  graphics 
routines  had  to  be  developed  for  the  system  used.  Also  the 
appropriate  variables  were  redimensioned  giving  consideration 
to  the  resolution  of  the  graphics  available  on  the 

10 


microcomputer  system  to  optimize  the  program  somewhat. 
Secondly,  although  it  is  possible  to  mix  graphics  and  text 
on  the  microcomputer  system,  it  is  not  done  in  a  straight  for¬ 
ward  manner  and  utilizes  extra  memory.  For  memory  considera¬ 
tions  and  convenience  it  was  decided  to  provide  sufficient 
information  for  interpretation  of  the  graphical  output 
below  each  graph.  It  is  felt  that  this  method  is  satisfac¬ 
tory  and  creates  little  inconvenience  to  the  user.  To 
enhance  the  interpretation  of  the  graphical  output  much 
time  and  effort  was  devoted  to  developing  plotting  routines 
to  display  the  data  in  relation  to  axes  labeled  with  tic-marks 
and  boundaries  of  integer  vice  fractional  values. 

A  major  obstacle  that  was  overcome  was  the  identification 
and  correction  of  a  memory  management  problem  unique  to  the 
microcomputer  system  used.  This  problem  affected  only  those 
programs  requiring  graphical  output.  The  essence  of  the 
problem  was  due  to  the  size  of  the  programs  and  the  number 
of  variables  used,  parts  of  the  program  and  stored  variables 
were  being  over  written  in  memory  by  the  graphics  routines. 
This  problem  was  finally  solved  by  making  appropriate  changes 
in  the  memory  management  scheme.  For  a  more  detailed  explana¬ 
tion  of  the  solution  see  Appendix  A. 


STATE  VARIABLE  PROGRAMS 


The  three  state  variable  programs  discussed  in  this  sec¬ 
tion  are  modified  versions  of  the  programs  of  the  same  names 
found  in  [Ref.  1].  These  programs  may  be  used  as  tools  for 
the  analysis  and  design  of  linear  control  systems  expressed 
in  the  following  state  variable  form: 


£(t)  =  Ax(t)  +  bu(t) 
u(t)  =  K[r(t )  -  kTx(t)] 
y(t)  =  cT£(t) 


The  Basic  Matrix  program  (BASMAT)  described  in  part  A 
of  this  section  is  used  to  compute  the  determinant,  inverse, 
characteristic  polynomial,  and  eigenvalues  of  the  square 
matrix  A.  It  is  also  used  to  determine  the  resolvent  matrix 
(sX-A)-1  and  state  transition  matrix  exp(At). 

The  Rational  Time  Response  program  (RTRESP),  described 
in  part  B,  is  used  to  determine  a  closed  form  expression 
for  the  time  response  of  a  system.  The  input  function  r(t) 
is  required  to  have  a  rational  time  response  and  no  repeated 
eigenvalues  are  allowed  in  the  combination  of  the  system 
and  input . 

The  Graphical  Time  Response  program  (RTRESP),  described 


in  part  C,  is  used  to  produce  a  graphical  display  of  the  time 
response  of  a  system  to  an  arbitrary  input. 


A.  BASMAT 


1.  Basic  Matrix  Program  ( BASMAT ) 

a.  Introduction 

When  done  by  hand,  matrix  manipulations  can  be 
quite  tedious  and  the  chances  of  an  error  being  made  are 
great.  In  the  study  of  linear  state  variable  analysis,  a 
computer  program  to  do  these  manipulations  is  almost 
essential.  It  can  do  the  necessary  manipulations  much  more 
quickly  and  accurately  and  allow  the  user  to  devote  his 
time  and  energy  to  design  and  analysis. 

b.  Description  of  Program 

BASMAT  [Ref.  1:  pp .  7,8]  takes  a  matrix  A, 
and  computes  the  determinant  of  A  (det  A),  the  inverse  of 
&  (A-1),  thr  characteristic  polynomial  (det (sI-£) ) ,  and 
eigenvalues  of  A  (Xf)as  well  as  the  resolvent  matrix 
0(s)  =  (sI-A)-1 
and  the  state  transition  matrix 
0(t)  *  exp (At) 

The  state  transition  matrix  is  expressed  as  matrix  coeffi¬ 
cients  times  the  natural  modes  exp(\t)  with  the  complex 
eigenvalues  expressed  as  damped  sine  and  cosine  terms.  The 
resolvent  matrix  is  written  as 
0(s)  =  adj (sI-A) /det (sI-A) 

and  the  matrix  numerator,  adj(sI-A),  is  output  as  matrix 
coefficients  of  powers  of  s  so  that  it  takes  the  form 
adj(sI-A)  =  F^  +  F£S  +  ...  +  F^s^-* 


The  BASMAT  program  interactively  accepts  input 
of  a  matrix  and  calls  subroutines  to  perform  the  appro¬ 
priate  calculations.  The  subroutines  used  are  CHREQ,  CHREQA , 
PROOT,  DET,  SIMEQ,  and  STMST.  These  subroutines  are  listed 
below  with  a  brief  description. 

CHREQ.  This  subroutine  is  used  to  determine  the 
characteristic  polynomial  det(sI-A),  and  the  resolvent 
matrix  (s^-^)-^  for  the  matrix  A.  The  Leverrier  algorithm 
is  used  to  compute  the  resolvent  matrix 
<£<s>  -  (Sl-A)*1 
in  the  form 

£<s>  = 

where  D(s)  is  the  characteristic  polynomial  det(sI-A). 

The  coefficients  of  the  characteristic  polynomial 
are  determined  by  subroutine  CHREQA. 

CHREQA.  This  subroutine  is  called  by  CHREQ  to 
determine  the  characteristic  polynomial,  det(sJ-A),  for  the 
A  matrix.  The  principal-minor  method  is  used. 

PROOT.  This  subroutine  uses  a  modified  Bair- 
stow  method  for  root  extraction  to  determine  the  roots  of  a 
polynomial  with  real  coefficients.  It  is  used  to  calculate 
the  eigenvalues  of  the  A  matrix. 

DET.  This  subroutine  computes  the  determinant 
of  a  matrix.  A  gauss  elimination  method  to  place  the  matrix 


N 


D(s) 
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in  upper  triangular  form  is  used.  It  is  called  by  CHREQA  to 
calculate  sub-determinants. 

SIMEQ.  This  subroutine  is  used  in  the  FORTRAN 
version  of  BASMAT  in  [Ref.  1:  pp.  7,8]  to  determine  the 
inverse  of  the  A  matrix.  In  the  BASIC  version  a  machine 
language  subroutine  [Ref.  2]  is  substituted  for  speed 
and  convenience. 

STMST.  This  subroutine  is  used  to  compute  the 
state  transition  matrix 
<^(t)  =  exp (At) 

for  a  matrix  A.  It  uses  the  Sylvester  Expansion  Theorem. 
Eigenvalues  of  the  A  matrix  must  be  provided.  This  routine 
can  not  handle  duplicate  eigenvalues. 

c.  Program  Translation  Problems 

See  section  II,  Translation  Considerations. 

2 .  BASMAT  User's  Guide 

This  program  occupies  about  7k  bytes  of  memory  and 
does  not  utilize  the  graphics  pages,  therefore  it  is  not 
necessary  to  relocate  the  disk  operating  system. 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

BASIC  MATRIX  PROGRAM 
PROBLEM  IDENTIFICATION  -> 


'6 


REMARKS : 


Label  the  problem.  Type  in  any  appropriate  combina¬ 
tion  of  letters,  numbers,  and/or  symbols,  excluding 
commas.  This  input  is  limited  to  255  characters. 

STEP  2: 

SCREEN  PROMPT: 

PLANT  ORDER?  -> 

REMARKS : 

Enter  the  number  that  corresponds  to  the  order  of 
the  A  matrix.  Maximum  order  allowed  is  ten. 

STEP  3: 

SCREEN  PROMPT: 

INPUT  PLANT  MATRIX. 

A(row, column )= 

REMARKS : 

Input  the  A  matrix  by  typing  in  each  element  as 
prompted.  The  program  will  ask  for  each  element 
beginning  with  the  first  row  and  going  from  left  to 
right . 

STEP  4: 

SCREEN  PROMPT: 

HARDCOPY?  (Y/N)  -> 

REMARKS : 

After  this  step  is  completed  the  program  will  output 
the  A  matrix  for  reference.  This  will  be  followed  by 
the  inverse  of  the  A  matrix,  the  determinant  of  the  A 


matrix,  the  matrix  coefficients  of  the  resolvent 
matrix  numerator,  the  characteristic  polynomial 
coefficients,  the  eigenvalues  of  the  plant  matrix, 
and  finally  the  elements  of  the  state  transition 
matrix.  See  Figure  1  for  a  sample  output. 

3.  BASMAT  Example 

This  example  will  use  the  system  matrix 


0  10 
A  =  0  0  1 

-2  -3  -3 


to  demonstrate  the  use  of  the  BASMAT  program.  It  will  refer 
to  the  steps  described  in  the  previous  section,  BASMAT 
User's  Guide. 

Step  1.  Enter  "EXAMPLE  <CR>"  where  <CR>  =  return. 

Step  2.  Enter  "3  <CR>" 

Step  3.  Enter  "0  <CR>  1  <CR>  0  <CR>  0  <CR>  0  <CR>  1  <CR>  -2 
<CR>  -3  <CR>  -3  <CR>" 

Step  4.  Enter  "Y  <CR>" 

The  resulting  output  is  shown  in  Figure  1.  These 
results  are  interpreted  as  follows: 

The  plant  matrix  is  shown  just  as  it  was  entered  with 
rows  horizontal  and  columns  vertical.  This  is  the  format 
used  for  all  matrix  output. 

The  inverse  of  the  plant  matrix  is  shown  followed  by 
the  scalar  determinant  of  the  A  matrix. 


The  matrix  coefficients  of  the  resolvent  matrix 
numerator  are  given  as  powers  of  s.  The  characteristic 
polynomial  is  listed  as  coefficients  of  powers  of  s. 

From  this  output ,  the  resolvent  matrix 
<^>(s)  =  adj(sI-A)/det(sI-£) 

may  be  written.  The  characteristic  polynomial,  det(sI-A), 

2  3 

is  2+3s+3s  +s  .  The  first  element  of  the  resolvent  matrix 

o 

numerator  is  3+3s+s  making  the  first  element  of  the 
resolvent  matrix 

\  3+3s+s 

<P//S)  - 2 3 

>l  2+3s+3s^+sJ 

The  real  and  imaginary  parts  of  the  eigenvalues  of 
the  plant  matrix  are  listed.  Finally,  the  elements  of  the 
state  transition  matrix  are  given.  The  first  element  of 
the  state  transition  matrix  can  be  written 

(t)  =  0.333e-2t+0.667e-0,5t  cos  0 . 866t+l .  15e"° ,5t  sin  0.866t 


BASIC  MATRIX  PROGRAM 
PROBLEM  IDENTIFICATION-EXAMPLE 

XXXXXXX*X**X*X*XXX*XX*XX**XXX**** 

THE  PLANT  MATRIX 

8  1  0 

0  0  1 

-2  -3  -3 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THE  INVERSE  OF  THE  PLANT  MATRIX 

-1.5  -1.5  -.5 

1  0  0 

-0  1  0 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THE  DETERMINANT  OF  THE  PLANT  MATRIX 
-2 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THE  MATRIX  COEFFICIENTS  OF  THE  RESOLVENT  MATRIX  NUMERATOR 


THE  MATRIX  COEFFICIENT  OF  SA2 


1  0  0 

0  1  0 

0  0  1 

THE  MATRIX  COEFFICIENT  OF  SA 1 

3  1  0 

0  3  1 

-2  -3  0 

THE  MATRIX  COEFFICIENT  OF  SA0 

3  3  1 

-2  0  0 

0-2  0 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX  Figure  1  BASMAT  Output 
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Figure  1  (Cont.) 

THE  CHARACTERISTIC  POLYNOMIAL  IN  ASCENDING  POWERS 
2  3  3  1 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
THE  EIGENVALUES  OF  THE  PLANT  MATRIX 
REAL  IMAGINARY 

-2  0 

-.5  -.866025484 

-.5  .366025404 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


THE  ELEMENTS  OF  THE  STATE  TRANSITION  MATRIX 


THE  MATRIX  COEFFICIENT  OF  EXP<-2)T 
.333333333  .333333333  .333333333 

— • 666666666  — . 666666666  — • 066606666 

1 . 33333333  1 . 33333333  1 . 33333333 

THE  MATRIX  COEFFICIENT  OF  EXP< -. 5) TSCOS< . 866025404) T 
.666666666  -.333333333  -.333333333 

•  666666666  1 • 0006600 7  • 666606060 

-1.33333333  -1.33333333  -.333333333 

THE  MATRIX  COEFFICIENT  OF  EXP< -. 5) TXSIN< . 8660 25404) T 
1.15470054  1.73205081  .577350263 

-1.15470054  -.577350268  0 

0  -1.15470054  -.577350268 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


Figure  1 


RTRESP 


1.  Rational  Time  Response  Program  (RTRESP) 

a.  Introduction 

Frequently  it  is  desirable  to  know  the  response 
of  a  system  as  a  function  of  time.  A  computer  program  can 
determine  this  quicker  and  more  accurately  than  by  hand. 

b.  Description  of  Program 

RTRESP  determines  the  time  response  in  closed 
form  of  the  closed  loop  system 

x(t )=Ax(t )+bu(t ) 
u(t )=K[r (t )-kTx(t ) ] 
y(t )=cTx(t ) 

due  to  any  initial  conditions  x(0)  and  input  r(t)  for  t>0. 
The  system  must  have  a  rational  Laplace  transform  R(s)  with 
a  pole-zero  excess  of  at  least  one.  [Ref.  1:  pp.  11,12] 

The  input  r(t)  is  treated  by  forming  a  m  -order 
dynamic  system  whose  initial  condition  response  is  equal  to 
r(t)  for  a  specific  set  of  initial  conditions.  This  system 
is  combined  with  the  original  system  and  then  the  complete 
response  in  closed  form  is  determined  from  the  subroutine 
STMST.  The  order  of  the  combined  system  must  be  ten  or 
less. 

Various  primary  and  utility  subroutines  are  used 
in  RTRESP.  The  primary  subroutines  are  listed  below  with  a 
brief  description: 


CHREQA.  This  subroutine  determines  the  charac¬ 
teristic  polynomial,  det(sI-A),  for  the  matrix  A  using  the 
principal-minor  method. 

PROOT.  This  subroutine  uses  a  modified  Bairstow 
method  for  root  extraction  to  determine  the  roots  of  a  poly¬ 
nomial  with  real  coefficients.  It  is  used  to  determine  the 
roots  of  the  numerator  and  denominator  of  R(s)  when  they  are 
entered  in  polynomial  form.  PROOT  is  also  used  to  determine 
the  eigenvalues  of  the  combined  system  matrix  (i.e.  the  roots 
of  the  characteristic  polynomial  determined  by  CHREQA). 

SEMBL.  This  subroutine  determines  the  coeffi¬ 
cients  of  a  polynomial  from  the  roots  of  the  polynomial. 

It  is  used  to  determine  the  polynomial  coefficients  of  the 
numerator  and  denominator  of  R(s)  when  they  are  entered 
in  the  factored  form.  This  subroutine  together  with  PROOT 
provide  the  feature  that  R(s)  may  be  entered  in  either 
factored  or  polynomial  form.  It  will  appear  in  the  output 
in  both  forms. 

STMST.  This  subroutine  computes  the  state  transi 

tion  matrix 
<j6(t)  =  exp(At) 

for  a  matrix  A.  It  uses  the  Sylvester  Expansion  Theorem. 
Eigenvalues  of  the  £  matrix  must  be  provided.  This  routine 
can  not  handle  duplicate  eigenvalues,  therefore  it  is 
necessary  that  the  combined  system  and  input  have  no 


repeated  eigenvalues.  STMST  is  used  to  determine  the  state 
transition  matrix  of  the  augmented  system, 
c.  Program  Translation  Problems 

See  section  II,  Translation  Considerations. 

2.  RTRESP  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

RATIONAL  TIME  RESPONSE  PROGRAM  (RTRESP) 

THIS  PROGRAM  DETERMINES  THE  TIME  RESPONSE  OF  A 
CLOSED-LOOP  SYSTEM  DUE  TO  SPECIFIED  INITIAL 
CONDITIONS  AND  INPUT.  SYSTEM  MUST  HAVE  A  RA¬ 
TIONAL  LAPLACE  TRANSFORM  WITH  A  POLE-ZERO  EXCESS 
OF  AT  LEAST  ONE. 

PROBLEM  IDENTIFICATION  -> 

REMARKS : 

Label  the  problem.  Type  in  any  appropriate  combina¬ 
tion  of  letters,  numbers,  and/or  symbols,  excluding 
commas  and  colons.  This  input  is  limited  to  255 
characters.  This  is  also  the  name  used  to  save  and 
retrieve  the  disk  file  containing  the  problem 
description.  (See  steps  2  and  3) 

STEP  2: 

SCREEN  PROMPT: 

WAS  THIS  PROBLEM  DESCRIPTION  PREVIOUSLY 
SAVED?  (Y/N)  -> 

REMARKS: 

This  step  provides  the  option  of  recalling  a  pre¬ 
viously  saved  problem  description  from  the  disk. 
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If  "Y"  is  typed,  the  program  looks  for  the  problem 
description  saved  under  the  label  typed  in  STEP  1. 
The  program  then  runs  to  completion  using  the 
retrieved  program  description.  If  an  "N"  is 
typed,  the  program  continues  with  STEP  3. 

STEP  3: 

SCREEN  PROMPT: 

DO  YOU  WANT  TO  SAVE  THIS  PROBLEM  DESCRIPTION? 
(Y/N)  -> 

REMARKS : 

A  positive  response  will  save  the  problem  descrip¬ 
tion,  that  will  be  entered  in  the  following  steps, 
to  the  disk  under  the  label  entered  in  STEP  1. 

This  problem  may  be  retrieved  during  a  later  session 
by  entering  the  proper  label  in  STEP  1  and  typing 
"Y"  in  STEP  2. 

STEP  4 : 

SCREEN  PROMPT: 

ORDER  OF  THE  SYSTEM?  -> 

REMARKS : 

The  total  order  of  the  system  and  the  input  must 
not  exceed  ten. 

STEP  5: 

SCREEN  PROMPT: 

INPUT  SYSTEM  (A)  MATRIX 


REMARKS : 


Enter  the  elements  of  the  A  matrix  as  prompted. 

The  format  is  A (row, column)  =  . 

STEP  6: 

SCREEN  PROMPT: 

THE  A  MATRIX 

(display  of  the  A  matrix) 

ANY  CHANGES?  (Y/N)  -> 

REMARKS : 

If  an  "NM  is  typed,  the  program  proceeds  to  STEP  9. 
To  correct  the  matrix  type  MY”  and  the  program  will 
proceed  to  STEP  7. 

STEP  7: 

SCREEN  PROMPT: 

TYPE  ROW, COLUMN  OF  THE  ELEMENT  TO  BE 
CORRECTED  -> 

STEP  8: 

SCREEN  PROMPT: 

A  (row, column)  = 

REMARKS : 

Enter  the  correct  value.  The  program  will  return 
to  STEP  6. 

STEP  9: 

SCREEN  PROMPT: 

INPUT  THE  CONTROL  (B)  VECTOR 


REMARKS : 


After  the  Jg  vector  is  entered,  a  correction  sequence 
similar  to  steps  6  through  8  is  encountered. 

STEP  10: 

SCREEN  PROMPT: 

INPUT  THE  OUTPUT  (C)  VECTOR 
C(l)  = 

REMARKS : 

After  the  C  vector  is  entered,  a  correction  sequence 
similar  to  steps  6  through  8  is  encountered. 

STEP  11: 

SCREEN  PROMPT: 

INPUT  THE  FEEDBACK  COEFFICIENTS 
FEEDBACK  COEFFICIENT  (1)  - 

REMARKS : 

After  the  feedback  coefficients  are  entered,  a 
correction  sequence  similar  to  steps  6  through  8 
is  encountered. 

STEP  12: 

SCREEN  PROMPT: 

INPUT  THE  GAIN  -> 

REMARKS : 


Enter  the  controller  gain 


STEP  13: 

SCREEN  PROMPT: 

INPUT  THE  INITIAL  CONDITIONS 
X0(1)  = 

REMARKS : 

After  the  initial  conditions  are  entered,  a 
correction  sequence  similar  to  steps  6  through 
8  is  encountered. 

STEP  14: 

SCREEN  PROMPT: 

ENTER  THE  INPUT  GAIN  -> 

REMARKS : 

Enter  the  gain  of  the  input  function,  R(s). 

STEP  15: 

SCREEN  PROMPT: 

'KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  NUMERATOR  POLYNOMIAL  -> 

REMARKS : 

Choose  the  preferred  method  of  entering  the  numerator 
polynomial  of  the  input  function,  R(f). 

STEP  16 A: 

SCREEN  PROMPT  (assumes  coefficient  form  chosen): 

INPUT  POLYNOMIAL  COEFFICIENTS 

INPUT  COEFFICIENTS  IN  ASCENDING  PWRS  OF  S 

INPUT  COEFFICIENT  OF  S^O  -> 

28 


H 


4 


A 


4 


1 


REMARKS : 


Enter  the  coefficient  along  with  its  algebraic 

sign  as  prompted.  The  program  assumes  that  the 

coefficient  of  the  highest  power  of  s  is  one. 

STEP  16B  (assumes  factored  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  FACTORS 

FOR  FACTOR  1 
REAL  PART  -> 

IMAGINARY  PART  -> 

REMARKS : 

Enter  the  factor  with  its  algebraic  sign.  Do  not 
enter  the  associated  root.  (e.g.  For  "(s-1)”  enter 

After  the  real  part  of  the  factor  is  entered, 
the  program  asks  for  the  imaginary  part,  allowing 
for  input  of  quadradic  factors.  If  a  non-zero 
imaginary  part  is  entered,  the  program  automatically 
enters  its  complex  conjugate. 

STEP  17: 

SCREEN  PROMPT: 

•KEY'  =  P  IOR  POLYNOMIAL  COEFFICIENT  FORM 
'KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  DENOMINATOR  POLYNOMIAL  -> 

REMARKS : 

The  denominator  is  entered  just  as  described  for 
the  numerator  in  steps  15  and  16. 


After  this  step  is  completed,  the  program  will 
output  all  input  for  reference  followed  by  the 
time  response  of  x(t)  and  the  system  output, 
y(t). 

3.  RTRESP  Example 

It  is  desired  to  know  the  time  response  of  the 
closed  loop  system 


x(t)  = 

"  0 

m 

1 

X  (t  )  + 

■*  «■ 

0 

-1 

-1_ 

_1 

u ( t )  =  5{r(t )  -  [1.0  0.5]  x ( t ) } 


y ( t )  =  [1  1]  x(t ) 


if  the  input  function  is  given  by 
R(s)  =  2.0  (1/s) 
and  the  initial  conditions  are 


x(t )  =  0 

Referring  to  the  steps  described  in  the  pervious  sec¬ 
tion,  RTRESP  User’s  Guide,  this  problem  would  be  entered  as 
follows : 

Step  1.  Enter  "EXAMPLE  <CR>"  where  <CR>  =  carriage 

return. 


Step  2  and  Step  3.  Enter  "N  <CR>." 

Step  4.  Enter  "2  <CR>." 

Step  5.  Enter  "0  <CR>  1  <CR>  -1  <CR>  -1  <CR>." 
Step  6.  Enter  "N  <CR>." 
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Step  9.  Enter  "0  <CR>  1  <CR> 


Enter  appropriate 


II 


response  for 

correction 

sequence . 

Step 

10. 

Enter 

"1  <CR>  1  <CR> . " 

Enter 

appropriate 

response  for 

correction 

sequence . 

Step 

11. 

Enter 

"1  <CR>  .5  <CR> . " 

Enter 

appropriate 

response  for 

correction 

sequence . 

Step 

12. 

Enter 

"5  <CR>." 

Step 

13. 

Enter 

"0  <CR>  0  <CR> . " 

Enter 

appropriate 

response  for 

correction 

sequence . 

Step 

14. 

Enter 

"2  <CR> . " 

Step 

15. 

Enter 

"P,0  <CR> . " 

Step 

16  A. 

Enter 

'  "1  <CR> . " 

Step 

17. 

Enter 

"F , 1  <CR>  0  <CR> . 

M 

The  program  will  now  run  to  completion  resulting  in  the  out¬ 
put  shown  in  Figure  2. 

The  output  begins  with  the  program  name  and  problem 
identification.  The  A,  B,  and  £  matrices  are  shown  for 
reference  followed  by  the  feedback  coefficients  and  the 
controller  gain.  Next,  the  initial  conditions  and  input 
gain  are  listed.  The  numerator  and  denominator  polynomials 
of  the  input  function  are  given  in  both  polynomial  coeffi¬ 
cient  form  and  factored  form  regardless  of  the  method  in 
which  they  were  entered.  The  polynomial  coefficient  form  is 
given  as  a  list  of  coefficients  from  left  to  right  in 
ascending  powers  of  s.  The  highest  power  of  s  is  always 
one.  The  polynomial  factored  form  appears  as  a  list  of  the 


real  and  imaginary  parts  of  each  root  of  the  polynomial. 

Notice  that  the  numerator  is  only  given  in  coefficient 
form  since  it  is  of  zero  order  and  no  roots  exist. 

The  time  response  of  the  state  x(t)  is  given  in  the 
form  of  vector  coefficients  of  the  various  natural  modes  of 
the  system  and  the  input  function.  From  this  output 
x(t)  is  seen  to  be 

X;L(t)  *  -1.67e”1,75t  cos  1 . 71t  -  1.70e"1,75t  sin  1.71t  +  1.67 
x2(t)  =  5.83e-1* 75t  sin  1.71t 

The  time  response  of  the  output  y(t)  is  given  as  scalor 
coefficients  of  the  same  natural  modes  as  x(t)  so  it  is  seen 
that 


=  -1 . 67e 


-1 . 75t 


1 


y(t ) 


cos  1.71t  +  4.13e 


-1 . 75t 


sin  1.71t  +  1.67 


RATIONAL  TIME  RESPONSE 

PROBLEM  IDENTIFICATION  ->  EXAMPLE 

THE  A  MATRIX 

0  1 

-1  -1 

THE  B  MATRIX 

0  1 

THE  C  MATRIX 

1  1 

FEEDBACK  COEFFICIENTS 
1  .5 

GAIN  =  5 

x*****xxxx**xxxs**x*r****x****x****rr****^*xs**sx 

INITIAL  CONDITIONS,  X<0> 

0  0 

RGAIN  =  2 


NUMERATOR  POLYNOMIAL  OF  R<S>  -  ASCENDING  PWRS  OF  S 
1 


Figure  2  RTRESP  Output 
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Figure  2  (Cont.) 


DENOMINATOR  POLYNOMIAL  OF  R<S>  -  ASCENDING  PWRS  OF  S 
9  1 

DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAGINARY  PART 

0  0 

***X**X**XXX**XX*****X**^**XX**********X*':^S*r*S 

THE  TIME  RESPONSE  OF  THE  STATE  X(T) 


THE  VECTOR  COEFFICIENT  OF  EXP< - 1 . 75) TXCOS< 1 . ? 1 3? 1365) T 
-1 .66666667  0 

THE  VECTOR  COEFFICIENT  OF  EXP< - 1 . 75) TSSIN< 1 . 7 13? 1365) T 
-1.70175824  5.83459946 

THE  VECTOR  COEFFICIENT  OF  EXP<0)T 
1.66466667  9 . 3 1322575E- 10 

XXX*XXXX*XXXXXXXXXXXXXXXXXX*XXXXXXr**XrA*r*Xr*XXXX 

THE  TIME  RESPONSE  OF  THE  OUTPUT  Y<T> 

THE  COEFFICIENT  OF  D<P< -  1 . 75) T2CQS<  1 . 7 139 1365) T 
-1 .66466467 

THE  COEFFICIENT  OF  EXP< -  1 . 75) T*SIN<  1 . 7 139 1365) T 
4.  13234143 

THE  COEFFICIENT  OF  EXP(0)T 


1 .66666447 


Figure  2 


GTRESP 


C. 

1.  Graphical  Time  Response  Program  (GTRESP) 

a.  Introduction 

Knowing  the  response  of  a  system  as  a  function  of 
time  can  be  very  helpful  in  the  study  and  analysis  of  that 
system.  The  graphical  display  of  this  response  can  give 
much  insight  into  a  system's  response  characteristics. 

b.  Description  of  Program 

GTRESP  [Ref.  1:  pp.  22-28]  determines  and 
graphically  displays  the  time  response  of  the  closed  loop 
system 

x(t)  =  Ax(t)  +  bu(t) 
u(t)  =  K  [r (t )  -  kT  x ( t )  ] 
y(t)  ■  cT  x(t) 

due  to  any  initial  conditions  x(0)  and  input  r(t). 

The  basic  purpose  of  this  program  is  similar 
to  that  of  the  RTRESP  program  described  earlier.  The 
difference  between  the  two  programs  is  that  GTRESP  will 
determine  the  time  response  for  arbitrary  input  functions 
which  may  not  have  rational  Laplace  transforms  but  RTRESP 
requires  a  rational  Laplace  transform.  Also  GTRESP  produces 
a  graphical  display  of  the  time  response  instead  of  a 
closed  form  expression  as  in  the  RTRESP  program. 

The  GTRESP  program  uses  a  fourth-order  Runge-Kutta 
numerical  integration  algorithm  to  calculate  the  time  response 


The  main  subroutines  used  are  CALCU,  RUNGE,  TRESP,  and  YDOT. 
Also  various  utility  and  plotting  subroutines  are  used. 

CALCU.  This  subroutine  is  called  by  TRESP  to 
determine  the  reference  input  r(t)  and  the  control  input 

u(t)=K  [r(t)  -  kTx(t)] 

The  reference  input  r(t)  must  be  defined  by  the  user  by 
inserting  the  appropriate  BASIC  coding  into  the  subroutine 
between  line  numbers  5010  and  5500.  The  reference  input  is 
represented  by  the  variable  R  and  the  control  input  by  the 
variable  U. 

RUNGE.  This  subroutine  is  called  by  TRESP  and 
contains  the  actual  fourth-order  Runge-Kutta  integration 
algorithm.  It  must  be  executed  four  times  for  each  inte¬ 
gration  step. 

TRESP.  This  subroutine  is  the  driving  subroutine 
which  calls  the  subroutines  CALCU,  RUNGE,  and  YDOT  along  with 
the  necessary  plotting  routines.  It  calculates  the  time 
response  of  the  closed  loop  linear  system  described  by  the 
input  parameters  and  plots  the  desired  variables. 

YDOT.  This  subroutine  is  called  by  TRESP  to 
compute  the  derivative 
x(t )=A£(t )+£u(t ) 

It  is  designed  to  handle  linear  systems  but  can  easily  be 
modified  to  handle  nonlinear  and  time  varying  systems  which 
would  give  GTRESP  a  nonlinear  and  time  varying  capability. 


36 


c.  Program  Translation  Problems 

See  section  II,  Translation  Considerations. 

2.  GTRESP  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

GRAPHICAL  TIME  RESPONSE  PROGRAM  (GTRESP) 

THIS  PROGRAM  DETERMINES  AND  GRAPHICALLY 
DISPLAYS  THE  TIME  RESPONSE  OF  A  CLOSED- 
LOOP  SYSTEM  DUE  TO  SPECIFIED  INITIAL 
CONDITIONS  AND  INPUT. 

PROBLEM  IDENTIFICATION  -> 

REMARKS : 

Label  the  problem.  Type  in  any  appropriate 
combination  of  letters,  numbers,  and/or  symbols, 
excluding  commas  and  colons.  This  input  is 
limited  to  255  characters.  This  is  also  the  name 
used  to  save  and  retrieve  the  disk  file  containing 
the  problem  description.  (See  steps  2  and  3) 

STEP  2: 

SCREEN  PROMPT: 

WAS  THIS  PROBLEM  DESCRIPTION  PREVIOUSLY 
SAVED?  (Y/N)  -> 

REMARKS : 

This  step  provides  the  option  of  recalling  a  pre¬ 
viously  saved  problem  description  from  the  disk. 

If  "Y"  is  typed  the  program  looks  for  the  problem 
description  saved  under  the  label  typed  in  STEP  1. 


The  program  then  runs  to  completion  using  the 
retrieved  program  description.  If  an  "N"  is 
typed  the  program  continues  with  STEP  3. 


STEP  3: 

SCREEN  PROMPT :  j 

DO  YOU  WANT  TO  SAVE  THIS  PROBLEM  DESCRIPTION? 

(Y/N)  -> 

REMARKS : 

I 

The  problem  description  will  be  entered  in  the 
following  steps. 

A  positive  response  will  save  the  problem  des- 

I 

cription  to  the  disk  under  the  label  entered  in 
STEP  1.  This  problem  may  be  retrieved  during  a 

later  session  by  entering  the  proper  label  in  STEP  j 

1  and  typing  "Y"  in  STEP  2. 

STEP  4: 

SCREEN  PROMPT: 

ORDER  OF  THE  SYSTEM?  -> 

REMARKS : 

The  total  order  of  the  system  and  the  input  must  not 
exceed  ten. 

STEP  5 : 

SCREEN  PROMPT:  ( 

INPUT  SYSTEM  (A)  MATRIX 


A(l,l)  = 


REMARKS 


Enter  the  appropriate  element  of  the  A  matrix.  The 
format  is  A(row, column )  =  . 

STEP  6: 

SCREEN  PROMPT: 

THE  A  MATRIX 

(display  of  the  A  matrix) 

ANY  CHANGES?  (Y/N)  -> 

REMARKS : 

If  an  "N"  is  typed,  the  program  proceeds  to  STEP  9. 
To  correct  the  matrix  type  "Y"  and  the  program  will 
proceed  to  STEP  7. 

STEP  7: 

SCREEN  PROMPT: 

TYPE  ROW, COLUMN  OF  THE  ELEMENT  TO  BE 
CORRECTED  -> 

STEP  8: 

SCREEN  PROMPT: 

A (row, column )  = 

REMARKS : 

Enter  the  correct  value.  The  program  will  return  to 
STEP  6. 

STEP  9: 

SCREEN  PROMPT: 

INPUT  THE  CONTROL  (B)  VECTOR 


REMARKS 


After  the  B  vector  is  entered  a  correction 
sequence  similar  to  steps  6  through  8  is 
encountered. 

STEP  10: 

SCREEN  PROMPT: 

INPUT  THE  OUTPUT  (C)  VECTOR 
C(l)  = 

REMARKS : 

After  the  C  vector  is  entered  a  correction  sequence 
similar  to  steps  6  through  8  is  encountered. 

STEP  11: 

SCREEN  PROMPT: 

INPUT  THE  FEEDBACK  COEFFICIENTS 
FEEDBACK  COEFFICIENT  (1)  = 

REMARKS : 

After  the  feedback  coefficients  are  entered  a  correc 
tion  sequence  similar  to  steps  6  through  8  is 
encountered. 

STEP  12: 

SCREEN  PROMPT: 

INPUT  THE  CONTROLLER  GAIN  -> 

REMARKS : 

Enter  the  controller  gain. 


STEP  13 


SCREEN  PROMPT: 


INPUT  THE  INITIAL  CONDITIONS 


XO(l)  = 


REMARKS : 


After  the  initial  conditions  are  entered  a  correc¬ 
tion  sequence  similar  to  steps  6  through  8  is 
encountered. 


STEP  14: 


SCREEN  PROMPT: 


ENTER  THE  FOLLOWING  PARAMETERS: 
NOTE. . (TF-TZ) / (DT*FR)  =<  100 

INITIAL  TIME  (TZ)  -> 


REMARKS : 


Enter  the  initial  time  of  the  time  interval  of 
interest.  Due  to  program  constraints  the  initial 
and  final  times  and  the  time  step  and  frequency  of 
output  must  be  chosen  so  that  they  satisfy  the 


relation. 


(TF-TZ)  /  (DT*FR)  =<  100 


STEP  15: 


SCREEN  PROMPT: 


FINAL  TIME  (TF)  -> 


REMARKS: 


Enter  the  final  time  of  the  time  interval  of 


interest 


STEP  16: 


SCREEN  PROMPT: 

TIME  STEP  (DT)  “> 


REMARKS : 

Enter  the  time  increment  for  each  step. 


STEP  17: 

SCREEN  PROMPT: 

FREQUENCY  OF  OUTPUT  (FR)  -> 


REMARKS : 

Enter  an  integer  n  and  the  program  will  print  out 
data  on  every  nth  time  step  iteration.  Ensure 
that  the  relation  expressed  in  the  remarks  of  STEP 
14  is  satisfied. 

STEP  18: 

SCREEN  PROMPT: 


YOU  MAY 

PLOT  UP  TO 

8  VARIABLES 

VS  TIME 

VARIABLE 

NUMBER 

VARIABLE 

NUMBER 

X1(T) 

1 

X8(T) 

8 

X2(T) 

2 

X9(T) 

9 

X3(T) 

3 

X10(T) 

10 

X4(T) 

4 

E(T) 

11 

X5(T) 

5 

U(T) 

12 

X6(T) 

6 

Y(T) 

13 

X7(T) 

7 

RCT) 

14 

HOW  MANY  VARIABLES  TO  PLOT?  MAX=8  “> 


REMARKS: 

You  may  plot  up  to  8  variables  vs  time.  These 
variables  will  be  plotted  on  the  same  plot  in  the 


output.  Enter  the  number  of  variables  you  want  to 
appear  on  the  plot. 


STEP  19: 

SCREEN  PROMPT: 

TYPE  THE  VARIABLE  NUMBER  <CR>  -> 

REMARKS : 

Type  the  number  associated  with  the  variable  of 
interest.  After  carriage  return  <CR>  is  typed,  the 
program  will  continue  to  prompt  for  another  variable 
number  until  the  number  of  variables  that  the  user 
indicated  in  STEP  18  has  been  entered. 

STEP  20: 

SCREEN  PROMPT: 

YOU  CHOSE  THE  FOLLOWING  VARIABLES: 

(list  of  variable  numbers  entered  in  STEP  19) 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES?  (Y/N)  -> 

REMARKS : 

To  make  a  change  in  the  variables  to  be  plotted 
type  a  "Y"  and  the  program  will  return  to  STEP 
18.  If  an  "N"  is  typed,  the  program  will  run  to 
completion. 

3.  GTRESP  Example 

It  is  desired  to  determine  and  plot  the  time  response 
of  the  error 

e(t)=r(t)-y(t) 

the  input  r(t),  and  the  state  variable  x0(t)  for  the  system 


x(t ) 


x(t)  + 


u(t ) 


1 


u(t)  =  5{r(t)  -  [1.0  0.5]  x(t)} 


y  (t )  *  [1  1]  x ( t ) 


The  time  interval  of  interest  is 
0<t  <4 


The  iteration  step  size  (DT)  is  chosen  to  be  0.01  and  printed 
output  is  desired  every  ten  steps  so  FREQ  equals  10.  The 
initial  conditions  are 
x(t )=0 

and  the  input  function  is 


r(t)  =  5.0  if  0<t<l 

r(t)  =  0  otherwise. 

In  order  to  define  this  input  function  the  following 


I 


BASIC  coding  is  inserted  between  line  numbers  5010  and  5500 

in  the  subroutine  CALC. 

5020  IF  TZ  >  1  GOTO  5040 
5030  R  =  5:  GOTO  5500 
5040  R  =  0 

5500  REM  END  OF  ROUTINE  DESCRIBING  R(T ) 


Referring  to  the  steps  described  in  section  2, 
GTRESP  User's  Guide,  this  problem  would  be  entered  as 
follows : 

Step  1.  Enter  "EXAMPLE  <CR>"  where  <CR>  =  carriage 


I 


return 


Step  2  and  Step  3.  Enter  "N  <CR>. 
Step  4.  Enter  "2  <CR>." 


Step  5.  Enter  "0  <CR>  1  <CR>  -1  <CR>  -1  <CR>." 

Step  6.  Enter  "N  <CR>." 

Step  9.  Enter  "0  <CR>  1  <CR>."  Enter  appropriate 
response  for  correction  sequence. 

Step  10.  Enter  "1  <CR>  1  <CR>."  Enter  appropriate 
response  for  correction  sequence. 

Step  11.  Enter  "1  <CR>  .5  <CR>."  Enter  appropriate 
response  for  correction  sequence. 

Step  12.  Enter  ”5  <CR>." 

Step  13.  Enter  "0  <CR>  0  <CR>."  Enter  appropriate 
response  for  correction  sequence. 

Step  14.  Enter  "0  <CR>." 

Step  15.  Enter  ”4  <CR>." 

Step  16.  Enter  ”.01  <CR>." 

Step  17.  Enter  "10  <CR>."  Ensure  that  the  constraint 
given  in  step  14 -of  GTRESP  User's  Guide  is  met.  That  con¬ 
straint  is  ( TF-TZ ) / ( DT*FR )  =  100 

In  this  case 
(4-0 ) / ( . 01*10 ) =40 
so  the  constraint  is  met . 

Step  18.  Enter  ”3  <CR>." 

Step  19.  Enter  "11  <CR>  14  <CR>  2  <CR>."  This  causes 
the  error  e(t),  the  input  r(t),  and  the  state  variable  x2(t), 
respectively,  to  be  plotted  as  desired  from  the  problem 
statement . 
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The  program  will  now  run  to  completion  resulting  in 


the  output  shown  in  Figure  3. 

The  output  begins  with  the  program  name  and  problem 
identification.  The  A,  B,  and  C  matrices  are  shown  for 
reference  followed  by  the  feedback  coefficients  and  the 
controller  gain.  Next  the  initial  conditions  and  time 
parameters  are  listed. 

The  second  page  of  output  lists  in  tabular  form  the 
value  of  time  t  and  the  corresponding  values  of  the  output 
y(t),  the  control  u(t),  and  all  of  the  state  variables. 

The  user  has  no  control  over  the  variables  output  in  this 
form. 

The  third  page  of  output  is  the  graph  itself.  Below 
the  graph  is  enough  information  to  properly  interpret  the 
plotted  data. 


GRAPHICAL  TIME  RESPONSE 

PROBLEM  IDENTIFICATION  ->  EXAMPLE 

************************************************* 

THE  A  MATRIX 

0  1 

-1  -1 

THE  B  MATRIX 
8  1 

THE  C  MATRIX 
1  1 

FEEDBACK  COEFFICIENTS 
1  .5 

GAIN  *  5 


INITIAL  CONDITIONS,  X  <  0  ) 
0  9 


INITIAL  TIME 


FINAL  TIME  =  4 


TIME  STEP  =  .91 

**************************** 


FREQUENCY  OF  OUTPUT 
*  *  *  *******  *  ********** 


Figure  3  GTRESP  Output 


Figure  3  (Cont.) 


GRAPHICAL  TIME  RESPONSE 
PROBLEM  IDENTIFICATION  ->  EXAMPLE 


T 

Y<T> 

IKT) 

X  1  <T) 

X2<T) 

0 

0 

25 

0 

8 

.1 

2.19944778 

19.2237187 

.111844751 

2.38838303 

.2 

3.34844475 

14.3947742 

.39362274 

3.454344 

.3 

5.0245394 

10.4744595 

.782794787 

4.24374241 

.4 

5.31339897 

7.39775741 

1.22749887 

4.5859889 

.5 

4.28493902 

5.0453327 

1.4389279 

4.59481112 

.599999998 

6.51052439 

3.37443244 

2.13898243 

4.37142176 

.499999996 

4.55142505 

2.22609087 

2.55813841 

3.99323444 

.799999994 

4.44012031 

1.51321108 

2.93459526 

3.52552585 

.899999991 

4.28025128 

1.14451516 

3.26194244 

3.31830842 

.999999989 

6.84703255 

1.93493374 

3.53319394 

2.58883361 

1.09999999 

3.42037346 

-18.1932886 

3.45493877 

-.0345453122 

1.19999998 

1.69848535 

-13.1375299 

3.5544866 

-1.85730124 

1.29999998 

.258794742 

-8.91302849 

3.38441444 

-3.84762187 

1.39999998 

-.748858826 

-5.48709854 

2.96349824 

-3.73255737 

1.49999998 

-1.4543045 

-2.79448401 

2.57289891 

-4.82728341 

1.59999998 

-1.84407253 

-.759574484 

2.16790317 

-4.83197569 

1.49999998 

-2.85893263 

.714871188 

1.77338415 

-3.83223678 

1.79999997 

-2.09253403 

1.71454348 

1.48590864 

-3.49844248 

1.89999997 

-2.01080576 

2.33453974 

1.87418985 

-3.88499541 

1.99999997 

-1.85192291 

2.45573142 

.78963026 

-2.64155317 

2.89999994 

-1.44447353 

2.7448857 

.547919244 

-2.19459277 

2.19999994 

-1.41908918 

2.47275259 

.34998815 

-1.76987733 

2.29999994 

-1.18722972 

2.48587785 

. 192873584 

-1.38010831 

2.39999994 

-.964037042 

2.22894883 

.8724495898 

-1.03648455 

2.49999996 

-.758194579 

1.93544872 

-.8140449082 

-.742129671 

2.59999995 

-.574947448 

1.4314455 

-.8774387299 

-.497316738 

2.69999995 

-.414853701 

1.33489444 

-.117184954 

-.299748744 

2.79999995 

-.284448821 

1.85847972 

-.139023864 

-.145425754 

2.89999995 

-.174814194 

.818704287 

-.147445518 

-.8293586731 

2.99999995 

-.8920420427 

.595121516 

-.145984544 

.853924521 

3.09999994 

-.0274993774 

.413224971 

-.137590611 

.189891233 

3.19999994 

.0198515208 

.244229347 

-.124743259 

.14379478 

3.29999994 

.0510344439 

.145923933 

-.109484217 

.168442841 

3.39999994 

.0710053207 

.8552049217 

-.8930872894 

.16409261 

3.49999993 

.081500273 

-.8115845112 

-.8748984485 

.158398742 

3.59999993 

.0847905445 

-.85793434 

-.0616148085 

.144487353 

3.49999993 

.882838023 

-.8877345728 

-.8477441939 

.138582217 

3.79999993 

.0772893216 

-.104315086 

-.8355433192 

.112352641 

3.89999993 

.8494859387 

-.118744482 

-.8251881781 

.8944741148 

3.99999992 

.0404847848 

-.109494415 

-.0166081408 

.3778949274 

Figure  3 
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ABSCISSA  ->  TIME  AXIS 
ORDINATE  ->  RESPONSE  MAGNITUDE 
TIC  MARKS  SHOW  INTERVALS  OF  UNITY 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  0  TO  A 
ORDINATE,  -3  TO  5 

SYSTEM  RESPONSE 

VARIABLE  SYMBOL 

X2  ♦ 

ERRORX 

INPUTT 


Figure  3 
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FRESP 


1.  Frequency  Response  Program  (FRESP) 

a.  Introduction 

The  response  of  a  system  as  a  function  of  fre¬ 
quency  is  a  very  important  characteristic  of  that  system. 

Some  common  ways  of  graphically  displaying  the  frequency 
response  of  a  system  include  the  use  of  amplitude  and  phase 
Bode  plots  and  Nyquist  diagrams.  These  are  valuable  tools 
for  system  analysis. 

b.  Description  of  Program 

The  FRESP  program  [Ref.  1:  pp.  105-113]  is 
used  to  determine  the  frequency  response  of  a  rational  trans¬ 
fer  function  G(s)  of  the  form 


G(s) 


„  A(s) 
K  B(s) 


where 

ACs)  =  +  a2s  +  ...  .  +  sM 

B(s)  *  b^  +  bgS  +  .  .  .  +  bjjS  +  s 

The  output  may  take  the  form  of  Bode  plots  or  a  Nyquist  dia¬ 

gram  or  both  in  addition  to  tabular  data. 

FRESP  gives  the  user  the  option  of  supplying 
discrete  frequency  values  or  allowing  the  program  to  linearly 
or  logarithmically  interpolate  frequency  values  between  two 
limit  values.  The  complex  number  G(jCJ)  can  be  computed  for 
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each  frequency  value.  G(j(jJ)  can  be  written  in  rectangular 
form  as 

G<jO»  =  R(  (jJ)  +  JX(CJ) 

where  R(£J)  and  X(GJ)  are  real  values.  The  magnitude  and 
phase  of  G(  j  (jJ )  can  then  be  written  as 

G(jCJ)  =  [R2(CJ)  +  X2((x))]i 
arg  G(j£J)  =  arctan  X(6J)/R(CJ) 

The  FRESP  program  uses  the  subroutines  PROOT, 

PVAL,  and  SEMBL.  The  subroutines  MAXI,  PHNOM,  GRAPH,  and 
SPLIT  that  appear  in  the  FORTRAN  version  do  not  appear  in  the 
BASIC  version.  They  are  incorporated  into  the  main  program 
and  into  various  plotting  subroutines  written  for  the  specific 
microcomputer  system  used. 

PROOT.  This  subroutine  uses  a  modified  Bairstow 
method  for  root  extraction  to  determine  the  roots  of  a 
polynomial  with  real  coefficients.  It  is  used  to  determine 
the  roots  of  the  numerator  and  denominator  of  G(s)  when 
they  are  entered  in  polynomial  form. 

SEMBL.  This  subroutine  determines  the  coeffi¬ 
cients  of  a  polynomial  from  the  roots  of  the  polynomial. 

It  is  used  to  determine  the  polynomial  coefficients  of  the 
numerator  and  denominator  of  G(s)  when  they  are  entered  in 
the  factored  form.  This  subroutine  together  with  PROOT 
provides  the  feature  that  G(s)  may  be  entered  in  either 
factored  or  coefficient  form  and  it  will  appear  in  the  output 
in  both  forms. 
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PVAL.  This  subroutine  is  used  to  evaluate  a 
polynomial  A(s)  with  real  coefficients  for 
s  =  PR  +  jPI 

PVAL  is  executed  twice  for  each  frequency  value  used,  once 
for  the  numerator  of  the  transfer  function  and  again  for 
the  denominator. 

c.  Program  Translation  Problems 

In  addition  to  the  programing  considerations  dis¬ 
cussed  in  section  II,  an  output  anomaly  was  traced  to  an 
apparent  oversight  in  the  subroutine  PVAL.  It  was  discovered 
that  a  zero  order  numeraxor  over  a  second  order  denomina¬ 
tor  with  a  free  s  such  as 
l/s(s+10) 
was  plotted  as 

1/ (s+10 ) 

The  reason  for  this  was  that  PVAL  treated  a  zero  order 

polynomial  and  a  first  order  polynomial  both  as  a  first  order 

polynomial  which  in  this  case  resulted  in  the  free  s  in 

the  denominator  being  cancelled.  The  reason  for  this  can 

be  seen  by  examining  the  portion  of  FORTRAN  code  from  PVAL 

in  [Ref.  1:  p.  164]  repeated  below: 

P=CMPLX(A(NN+1) ,0.  ) 

DO  100  J=1,NN 
100  P=P*S+A(NN+1-J') 

VR=REAL(P) 

In  this  portion  of  code  NN  represents  the  order 
of  the  polynomial  being  evaluated  and  the  intention  is  that 
the  DO  loop  be  executed  a  number  of  times  equal  to  the  order 


of  the  polynomial.  But  because  of  the  nature  of  the  DO  loop, 

if  the  polynomial  is  of  zero  order  (i.e.  NN  equals  zero) 

the  loop  will  be  executed  once  just  as  if  the  polynomial 

was  of  first  order  [Ref.  3].  One  way  to  correct  this 

problem  is  to  modify  the  section  of  PVAL  as  shown  here: 

P=CMPLX( A(NN+1 ) ,0. ) 

IF  (NN)  110,110,90 
90  DO  100  J=1 ,NN 
100  P=P*S+A(NN+1-J) 

110  VR=REAL(P) 

The  essence  of  the  modification  is  to  skip  the  DO  loop  if 
the  polynomial  has  order  zero.  Translating  this  modifica¬ 
tion  into  the  BASIC  version  of  FRESP  solved  the  problem. 

2 .  FRESP  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

FREQUENCY  RESPONSE  PROGRAM  (FRESP) 

THIS  PROGRAM  OBTAINS  AND  PLOTS  THE  FREQUENCY 
OF  A  RATIONAL  TRANSFER  FUNCTION  OVER  A  SPECIFIED 
RANGE  OF  FREQUENCIES.  BOTH  RECTANGULAR  BODE 
PLOTS  AS  WELL  AS  A  POLAR  NYQUIST  PLOT  CAN  BE 
OBTAINED. 

PROBLEM  IDENTIFICATION  “> 

REMARKS : 

Label  the  problem.  Type  in  any  appropriate  combina¬ 
tion  of  letters,  numbers,  and/or  symbols,  excluding 
commas  and  colons.  This  input  is  limited  to  255 
characters. 
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STEP  2 


SCREEN  PROMPT: 

INPUT  THE  GAIN  -> 

REMARKS : 

Enter  the  transfer  function  gain. 

STEP  3: 

SCREEN  PROMPT: 

•KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY*  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  NUMERATOR  ,A(S)  -> 

REMARKS : 

Choose  the  preferred  method  of  entering  the  numera¬ 
tor  polynomial  of  the  transfer  function,  G(s). 

STEP  4A  (assumes  coefficient  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  COEFFICIENTS 

INPUT  COEFFICIENTS  IN  ASCENDING  PWRS  OF  S 

INPUT  COEFFICIENT  OF  S~0  -> 

REMARKS : 

Enter  the  coefficient  along  with  its  algebraic  sign 
as  prompted.  The  program  assumes  that  the  coeffi¬ 
cient  of  the  highest  power  of  s  is  one. 

STEP  4B  (assumes  factored  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  FACTORS 

FOR  FACTOR  1 
REAL  PART  -> 


IMAGINARY  PART  -> 


REMARKS : 

Enter  the  factor  with  its  algebraic  sign.  Do  not 
enter  the  associated  root.  (e.g.  For  "(s-l)" 
enter  "-1".)  After  the  real  part  of  the  factor  is 
entered,  the  program  asks  for  the  imaginary  part 
allowing  for  input  of  quadradic  factors.  If  a 
non-zero  imaginary  part  is  entered,  the  program 
automatically  enters  its  complex  conjugate. 

STEP  5: 

SCREEN  PROMPT: 

■KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  *  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  DENOMINATOR,  B(S)  -> 

REMARKS : 

The  denominator  is  entered  just  as  described  for 
the  numerator  in  steps  3  and  4. 

STEP  6: 

SCREEN  PROMPT: 

MINIMUM  FREQUENCY  VALUE  -> 

REMARKS : 

Enter  the  minimum  frequency  value  in  radians  ,er 
second  of  the  frequency  range  of  interest. 

STEP  7: 

SCREEN  PROMPT: 

MAXIMUM  FREQUENCY  VALUE  -> 


REMARKS : 


Enter  the  maximum  frequency  value  in  radians  per 
second  of  the  frequency  range  of  interest. 

STEP  8: 

SCREEN  PROMPT: 

NUMBER  OF  FREQUENCY  VALUES  TO  BE  USED  -> 

REMARKS : 

Enter  an  integer  from  1  to  200.  Since  the  horizon¬ 
tal  resolution  capability  of  the  graphics  on  the 
microcomputer  system  used  is  191  pixels,  any  number 
of  values  greater  than  191  adds  no  more  detail  to  the 
Bode  plots.  For  the  Nyquist  diagram  the  greater  the 
number  of  values  the  greater  the  detail  due  to  the 
nature  of  the  plot.  If  a  list  of  discrete  fre¬ 
quency  values  will  be  supplied  (i.e.  option  1  in 
step  9)  enter  the  number  of  frequencies  in  that 
list . 

STEP  9: 

SCREEN  PROMPT: 

0  -  LOGARITHMIC  INTERPOLATION 

1  =  DISCRETE  VALUES  SUPPLIED 

2  =  LINEAR  INTERPOLATION 

CHOOSE  ONE  -> 

REMARKS : 

Make  a  choice  by  entering  the  associated  number.  This 
step  allows  the  user  to  choose  the  method  used  by  the 
program  to  obtain  discrete  frequency  values  used  to 
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generate  data  points.  LOGARITHMIC  INTERPOLATION  pro¬ 
duces  data  points  that  appear  equally  spaced  on  a 
logarithmic  scale  between  the  minimum  and  maximum 
frequency  values  entered  in  steps  6  and  7.  With  this 
choice  a  Bode  and/or  Nyquist  plot  may  be  obtained. 
LINEAR  INTERPOLATION  produces  data  points  that  appear 
equally  spaces  on  a  linear  scale.  Only  a  Nyquist 
plot  may  be  obtained  with  this  choice.  DISCRETE  VALUES 
SUPPLIED  allows  the  user  to  supply  a  list  of 
discrete  frequency  values  of  interest.  No  plots 
may  be  obtained,  however,  with  this  option.  Below 
is  a  table  summarizing  the  available  graphical  output 
for  each  choice. 


CHOICE 

BODE 

NYQUIST 

0 

YES 

YES 

1 

NO 

NO 

2 

NO 

YES 

STEP  10  (assumes  choice  1  in  step  9): 

SCREEN  PROMPT: 

TYPE  IN  THE  DISCRETE  FREQUENCY  VALUES, 

(number  entered  in  step  8)  VALUES  NEEDED 
FREQ  ( 1 ) ? 

REMARKS : 

Enter  the  list  of  frequency  values  as  prompted.  When 
the  last  value  has  been  entered  the  program  will  run 
to  completion  outputting  data  similar  to  Figure  4  with 
the  exception  of  the  plots. 


STEP  11  (assumes  choice  0  in  step  9): 

SCREEN  PROMPT: 

BODE  PLOT?  (Y/N)  -> 

REMARKS : 

The  program  outputs  tabular  data  regardless  of  the 
choice  made  here. 

STEP  12  (assumes  choice  0  or  2  in  step  9): 

SCREEN  PROMPT: 

NYQUIST  PLOT?  (Y/N)  -> 

REMARKS : 

The  program  outputs  tabular  data  regardless  of  the 
choice  made  here. 

After  this  step,  the  program  runs  to  completion. 

See  sample  output  in  Figure  4. 

3.  FRESP  Example 

It  is  desired  to  know  the  frequency  response  of  the 
transfer  function 

8(0.5  +  s) 

G(s)  =  - 

4  +  6s  +  3s2  +  s3 


The  frequency  range  of  interest  is  CJ  =  0.1  rad/sec  to  100 
rad/sec.  A  Bode  plot  of  the  magnitude  and  phase  is  desired 
but  a  Nyquist  plot  is  not.  The  plots  should  be  generated 
from  100  values  of  frequency  logarithmically  spaced  from 
0.1  to  100. 


Referring  to  the  steps  described  in  the  previous 


section 

,  FRESP  User's  Guide,  this  problem  would  be  entered 

as  follows: 

STEP 

1. 

Enter 

"EXAMPLE  CR  "  where  <CR>  =  carriage 

return. 

STEP 

2. 

Enter 

"8  <CR> . " 

STEP 

3. 

Enter 

"F , 1  <CR> . " 

STEP 

4B . 

Enter 

".5  <CR>  0  <CR> . " 

STEP 

5. 

Enter 

"P , 3  <CR>  4  <CR>  6  <CR>  3  <CR>." 

STEP 

6. 

Enter 

".1  <CR> . " 

STEP 

7. 

Enter 

"100  <CR> . " 

STEP 

8. 

Enter 

"100  <CR> . " 

STEP 

9. 

Enter 

"0  <CR> . " 

STEP 

11. 

Enter 

"Y  <CR> . " 

STEP 

12. 

Enter 

"N  <CR> . " 

The  program  will  then  continue  to  completion  pro¬ 
ducing  the  output  seen  in  Figure  4. 

The  output  begins  with  the  program  name  and  problem 
identification.  The  transfer  function  gain  is  given  followed 
by  the  transfer  function  numerator  and  denominator  each 
listed  as  coefficients  in  ascending  powers  of  s  and  then  as 
the  real  and  imaginary  parts  of  the  roots. 

The  second  page  of  output  is  tabular  data.  It 
consists  of  the  radian  frequency  and  the  transfer  function's 
corresponding  real  and  imaginary  parts,  magnitude,  and  phase 
in  radians  and  in  degrees.  Although  100  frequency  values 


were  generated  and  their  associated  data  points  plotted, 
tabular  data  appears  only  for  every  other  data  point.  To 
reduce  unnecessary  output  and  for  formatting  purposes  the 
following  scheme  is  used  for  tabular  data. 

NUMBER  OF  FREQUENCIES  TABULAR  DATA 

REQUIRED  FROM  STEP  8  PRINTED  FOR 

0  to  50  every  freq. 

51  to  100  every  other  freq. 

101  to  150  every  3rd  freq. 

151  to  200  every  4th  freq. 

The  next  two  pages  of  output  are  the  Bode  plots  for 
amplitude  and  phase.  The  plots  are  headed  with  the  problem 
identification  and  type  of  plot  and  below  each  plot  is 
sufficient  information  to  allow  proper  interpretation  of  the 
plotted  data.  Note  that  the  phase  angles  are  normalized  to 
always  remain  between  -180  and  +180  degrees. 

This  is  the  end  of  the  output  generated  from  the 
example  as  it  was  input  since  a  Nyquist  plot  was  not  desired 
A  Nyquist  plot  is  included  as  the  last  page  of  Figure  4  for 
the  sake  of  completeness. 


FREQUENCY  RESPONSE 

PROBLEM  IDENTIFICATION  -  EXAMPLE 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*** 

GAIN=8 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

.5 

1 


NUMERATOR  ROOTS  ARE 
REAL  PART 


IMAGINARY  PART 


-.5 


d 


DENOMINATOR  COEFFICIENTS 

4 

6 

3 

1 

DENOMINATOR  ROOTS  ARE 
REAL  PART 

-1 

-1 

-1 


IN  ASCENDING  POWERS  OF  S 


IMAGINARY  PART 

-1.73205081 

1.73205081 

0 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


Figure  4  FRESP  Output 


Figure  4 

(Cont.  )  PROBLEM  IDENTIFICATION 

-  EXAMPLE 

RADIAN  FREQ. 

REAL  PART 

IMAGINARY  PART 

MAGNITUDE 

PHASE (RAD) 

PHASE  (CEu) 

.187226722 

1.81785827 

.8516283768 

1.81836782 

.8587189222 

2.98593123 

.123284674 

1.82247187 

.8586342813 

1.82415889 

.8572828479 

3.<.328666 

.141747416 

1.82956848 

.8663192682 

1.83178223 

.8643257561 

3.68559566 

.162978884 

1.83885884 

.874688284 

1.84152594 

.8716873348 

4.1873832 

.187381743 

1.85894983 

.8832923195 

1.85424531 

.8798898887 

4.53146382 

.215443469 

1.86665859 

.8928537766 

1.87861542 

.8868884184 

4.93258435 

.247787636 

1.88698827 

.188388178 

1.89152785 

.8928269474 

5.27275758 

.284883588 

1.11285844 

.187136157 

1.11799565 

.8959768558 

5.4998249 

.327454917 

1.14575167 

.111131384 

1.1511286 

.896691838 

5.54883621 

.376493582 

1.18695955 

.118216475 

1.19286571 

.8925986297 

5.3858542 

.43287613 

1.23774461 

.181419568 

1.24189277 

.8817563658 

4.68429638 

.497782359 

1.29983767 

.8886188318 

1.38153689 

.9619889349 

3.55124725 

.572236769 

1.37188815 

.8422676713 

1.37165954 

.3388198662 

1.76534389 

.657933228 

1.45241893 

-.8288852383 

1.45256188 

-.8143787124 

-.822839829 

.756463332 

1.53955782 

-.118868565 

1.54487711 

-.8765349743 

-4.38513258 

.869749888 

1.62457894 

-.268681164 

1.6453479 

-.159856477 

-9.11326889 

1.88888881 

1.69238769 

-.461538466 

1.75411684 

-.266252852 

-15.2551243 

1.149757 

1.71476725 

-.732787826 

1.36478888 

-.483858578 

-23.1389419 

1.32194116 

1.64366872 

-1.87486971 

1.96392254 

-.579143836 

-33.1825894 

1.51991189 

1.4183691 

-1.45816618 

2.82289964 

-.799389724 

-45.7978981 

1.74752841 

.965794858 

-1.7451883 

1.99459581 

-1.36534155 

-61.3395965 

2.88923382 

.388589559 

-1.7966967 

1.83654746 

-1.36289719 

-78.8424484 

2.31812972 

-.134536879 

-1.5559564 

1.56176192 

-1.65784786 

-94.9418367 

2.65688781 

-.423812889 

-1.16877489 

1.24324255 

-1.91866312 

-189.931338 

3.85385554 

-.58737737 

-.883954387 

.958678461 

-2.13377179 

-122.256162 

3.51119177 

-.478886616 

-.52956775 

.713394931 

-2.38586561 

-132.870578 

4.8378173 

-.486613557 

-.342879268 

.531884177 

-2.44182781 

-139.868595 

4.64158888 

-.32894126 

-.228992452 

.396282748 

-2.55882432 

-146.185683 

5.33669929 

-.259188135 

-.142546166 

.295723333 

-2.63862774 

-151.182287 

6.13598734 

-.281847626 

-.8921935898 

.221178189 

-2.71163856 

-155.3655 

7.85488239 

-.154626892 

-.85981862 

.165791393 

-2.77251189 

-158.353287 

3.1113884 

-.118278988 

-.8389127642 

.124515475 

-2.32375481 

-161.789291 

9.32683358 

-.898161962 

-.8253788175 

.8936654855 

-2.86722885 

-164.279667 

18.7226724 

-.8685723775 

-.8165836577 

.8785491932 

-2.98438789 

-166.484598 

12.3284676 

-.8528722615 

-.8188538543 

.8531914145 

-2.9368967 

-168.226009 

14.1747418 

-.8395884317 

-7.1 1247 156E-83 

.848135662 

-2.96344186 

-169.792726 

16.2975886 

-.829941485 

-4.66519619E-83 

.8383826696 

-2.98782466 

-171.143967 

18.7381745 

-.8226835836 

-3.8622238E-83 

.8228892671 

-3.8874861 

-172.311739 

21.5443472 

-.8171782888 

-2.81116278E-83 

.8172956174 

-3.8258474 

-173.322511 

24.7787639 

-.81388548 96 

-1.32142716E-03 

.8138724492 

-3.84833486 

-174.198418 

28.4883591 

-9.84425637E-83 

-8.68521676E-84 

9. 8824953 IE-83 

-3.85359432 

-174.958129 

32.7454921 

-7, 458269 18E-83 

-5.78986284E-84 

7.47211725E-83 

-3.86518255 

-175.617583 

37.6493587 

-5.63781716E-83 

-3.75458885E-84 

5.65838488E-83 

-3.87589595 

-176.198883 

43.2876135 

-4.26591987E-83 

-2.46911 164E-84 

4.2738587E-83 

-3.88377726 

-176.687485 

49.7782365 

-3.22764824E-83 

-1.62396323E-84 

3.23173I88E-83 

-3.89132896 

-177.119788 

57.2236776 

-2.44196898E-83 

-1.86818684E-84 

2.444296 15E-83 

-3.89787757 

-177.495374 

65.7933236 

-1. 84746 146E-83 

-7.82659939E-85 

1.84879722E-83 

-3.18357722 

-177.321939 

75.6463342 

-1.3976554 IE-83 

-4.62235779E-85 

1.39841956E-83 

-3.18853252 

-178.185858 

86.9749819 

-1.85734283E-83 

-3.84886657E-85 

1. 8577888  IE-83 

-3.11284111 

-173.352722 

188.888882 

-7.99879835E-84 

-2.88851985E-85 

8.88129964E-84 

-3.11658766 

-173.567383 

Figure  4 
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PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  BODE  PLOT  < AMPLITUDE)  XX 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 

MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  *  100  RADIANS/SEC 
AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-80  DECIBELS 


Figure  4 


Figure  4  (Cont.) 


PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  BODE  PLOT  (PHASE)  XX 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  98  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  180  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  *  90  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  =»  -188  DEGREES 


Figure  4 


ABSCISSA  ->  REAL  PART  OF  G<JW> 
ORDINATE  ->  IMAGINARY  PART  OF  G<JW> 
TIC  MARKS  SHOW  INTERVALS  OF  UNITY 
AXES  CROSS  AT  ORIGIN 


Figure  4 


RTLOC 


1.  Root  Locus  Program  (RTLOC) 

a.  Introduction 

When  analyzing  a  system  it  is  helpful  to  know  the 
location  of  the  closed  loop  poles  in  the  s-plane.  These 
closed  loop  poles  are  the  roots  of  the  characteristic 
equation  and  determine  the  basic  characteristics  of  the 
transient  response  of  a  closed  loop  system.  A  root  locus 
plot  is  a  plot  of  the  roots  of  the  characteristic  equation, 
usually  as  a  function  of  the  gain  of  the  transfer  function, 
and  therefore  it  is  a  valuable  tool  for  system  analysis. 

b.  Description  of  Program 

The  RTLOC  program  [Ref.  1:  pp.  114-121]  calcu¬ 
lates  and  plots  the  roots  of  the  equation 
1  +  KG(s )  =  0 

as  a  function  of  K.  G(s)  is  assumed  to  be  a  rational  function 
of  the  form 


G(s) 


N(s ) 
D(s) 


and  the  root  locus  becomes  the  locus  of  roots  of  DR(s)  as  K 
varies  where 

Dk(s)  =  DCs)  +  KN(s) 

RTLOC  uses  an  algebraic  plus  linear  progression 

scheme  to  vary  K  to  give  reasonable  spacing  of  the  roots. 

D^Cs)  is  obtained  for  each  value  of  K  and  the  subroutine  PROOT 
K. 


is  used  to  calculate  its  roots. 


The  scheme  used  to  calculate  values  of  K  assumes 


that  K  takes  on  only  positive  values.  If  K  is  to  range 
through  negative  value,  the  value  of  smaller  magnitude  (less 
negative)  a-ist  be  used  as  the  minimum  value.  The  routine 
starts  at  the  maximum  (less  negative)  value  and  becomes 
increasingly  negative  until  the  lower  limit  is  reached. 

If  both  positive  and  negative  values  are  desired  then  two 
separate  runs  must  be  made  with  only  positive  and  only 
negative  values. 

The  RTLOC  program  has  a  feature  that  allows  the 
user  to  specify  a  range  of  s  and  Cl)  values  around  an  area 
of  the  root  locus  plot  that  is  of  interest.  That  area  is 
then  enlarged  and  more  closely  spaced  values  of  K  are  generated 
giving  more  detail  to  the  plot. 

The  major  subroutines  used  are  PROOT  and  SEMBL. 

The  subroutine  SPLIT  used  in  the  FORTRAN  version  is  replaced 
by  additional  coding  in  the  main  program  and  by  various 
plotting  subroutines. 

PROOT.  This  subroutine  uses  a  modified  Bairstow 
method  for  root  extraction  to  determine  the  roots  of  a 
polynomial  with  real  coefficients.  It  is  used  to  determine 
the  roots  of  the  numerator  and  denominator  of  G(s)  when  it 
is  entered  in  polynomial  form.  It  is  also  used  to  determine 
the  roots  of  D^(s)  for  each  value  of  K  generated. 

SEMBL.  This  subroutine  determines  the  coeffi¬ 
cients  of  a  polynomial  from  the  roots  of  the  polynomial.  It 
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is  used  to  determine  the  polynomial  coefficients  of  the 
numerator  and  denominator  of  G(s)  when  they  are  entered  in 
the  factored  form.  This  subroutine  together  with  PROOT 
provides  the  feature  that  G(s)  may  be  entered  in  either 
factored  or  coefficient  form.  It  appears  in  the  output  in 
both  forms. 

c.  Program  Translation  Problems 

In  addition  to  the  programing  considerations 
discussed  in  section  II,  two  problems  were  encountered. 

When  a  range  of  negative  gains  is  entered  the  program  cal¬ 
culates  only  one  value  of  K  due  to  the  logic  of  line  229  in 
the  FORTRAN  code.  This  works  correctly  only  for  positive 
ranges  of  gain.  (Note  that  K  is  represented  by  the  variable 
G  in  the  computer  program. )  In  the  BASIC  version  coding 
was  added  to  test  the  sign  of  the  range  of  gains  and  modify 
the  logic  to  correctly  handle  the  negative  case. 

The  second  problem  was  encountered  when  using 
gains  of  very  small  magnitude  necessary  in  using  this 
program  for  w' -plane  analysis.  The  schemes  used  to  generate 
values  of  K  for  a  reasonable  spacing  of  the  roots  when  plotted 
worked  fine  for  the  magnitudes  usually  encountered  in  s-plane 
analysis  but  was  not  flexible  enough  to  handle  gains  of  much 
smaller  magnitudes  (e.g.  0  to  IE-3).  To  add  this  necessary 
flexibility  the  following  FORTRAN  lines  used  to  generate 
values  of  gain  G, 


227  G=1 . 15* (G+SIGNG*0 . 05 ) 

228  G=1 . 04* (G+SIGNG*0 . 02 ) 


were  modified  by  replacing  the  constant  values  0.05  and  0.02 

with  the  variables  D1  and  D2 ,  respectively,  where 

D1=ABS(GMIN-GMAX) /700 
D2=ABS(GMIN-GMAX) /1500 

so  that  the  gain  increment  was  a  function  of  the  range  of 
gains  of  interest . 

2 .  RTLOC  User's  Guide 

Push  "return"  after  each  input. 

STEP  1: 

SCREEN  PROMPT: 

ROOT  LOCUS  PROGRAM  (RTLOC) 

THIS  PROGRAM  PLOTS  THE  ROOT  LOCUS  OF  A 
DESCRIBED  SYSTEM. 

PROBLEM  IDENTIFICATION  -> 

REMARKS : 

Label  the  problem.  Type  in  any  appropriate  combina¬ 
tion  of  letters,  numbers,  and/or  symbols,  excluding 
commas  and  colons.  This  input  is  limited  to  255 
characters. 

STEP  2: 

SCREEN  PROMPT: 

INPUT  THE  RANGE  OF  GAINS  (MIN, MAX)  -> 

REMARKS : 

Enter  the  range  of  transfer  function  gains  of  in¬ 
terest.  Note  that  the  minimum  gain  is  the  gain  of 
lowest  absolute  value. 


STEP  3: 


SCREEN  PROMPT: 


'KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
'KEY*  =  F  FOR  POLYNOMIAL  FACTORED  FORM 


INPUT  'KEY',  ORDER  FOR  NUMERATOR,  N(S)  -> 

REMARKS : 

Choose  the  preferred  method  of  entering  the 
numerator  polynomial  of  the  transfer  function 
G(s). 

STEP  4A  (assumes  coefficient  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  COEFFICIENTS 

INPUT  COEFFICIENTS  IN  ASCENDING  PWRS  OF  s  OR 

w' 


INPUT  COEFFICIENT  OF  S~0  -> 

REMARKS : 

Enter  the  coefficient  along  with  its  algebraic  sign 
as  prompted.  The  program  assumes  that  the 
coefficient  of  the  highest  power  of  s  is  one. 

STEP  4B  (assumes  factored  form  chosen): 

SCREEN  PROMPT: 

INPUT  POLYNOMIAL  FACTORS 


FOR  FACTOR  1 
REAL  PART  -> 


IMAGINARY  PART  -> 


REMARKS : 


Enter  the  factor  with  its  algebraic  sign.  Do  not 
enter  the  associated  root.  (e.g.  For  "(s-1)"  enter 
"-1" .)  After  the  real  part  of  the  factor  is 
entered  the  program  asks  for  the  imaginary  part 
allowing  for  input  of  quadradic  factors.  If  a 
non-zero  imaginary  part  is  entered  the  program 
automatically  enters  its  complex  conjugate. 

STEP  5: 

SCREEN  PROMPT: 

'KEY'  =  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
•KEY'  =  F  FOR  POLYNOMIAL  FACTORED  FORM 

INPUT  'KEY',  ORDER  FOR  DENOMINATOR,  DCS)  “> 

REMARKS : 

The  denominator  is  entered  just  as  described  for  the 
numerator  in  steps  3  and  4. 

STEP  6: 

SCREEN  PROMPT: 

WOULD  YOU  LIKE  TO  LOOK  AT  ONLY  A  PART  OF  THE 
ROOT  LOCUS?  (Y/N)  -> 

REMARKS : 

Enter  an  "N"  for  the  option  of  viewing  the  root  locus 
for  the  entire  range  of  gains  entered  in  STEP  2. 

Enter  a  "Y"  for  the  option  of  viewing  only  a  portion 
of  the  root  locus  defined  by  minimum  and  maximum 
values  of  sigma  and  OJ  of  interest  to  the  user. 
Usually,  it  is  best  to  view  the  entire  root  locus 
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first  then  decide  if  any  section  needs  to  be 
enlarged  to  reveal  more  detail  and  rerun  RTLOC 
to  view  only  that  portion. 

STEP  7  (assumes  "Y"  entered  in  step  6): 

SCREEN  PROMPT: 

ENTER  SIGMA  MIN,  SIGMA  MAX  -> 

REMARKS : 

Enter  the  minimum  and  maximum  values  of  sigma  (i.e. 
the  real  axis)  of  the  portion  of  the  root  locus 
plot  of  interest. 

STEP  8  (assumes  "Y"  entered  in  step  6): 

SCREEN  PROMPT: 

ENTER  OMEGA  MIN,  OMEGA  MAX  -> 

REMARKS : 

Enter  the  minimum  and  maximum  values  of  omega  (i.e. 
the  imaginary  axis)  of  the  portion  of  the  root 
locus  plot  of  interest. 

STEP  9: 

SCREEN  PROMPT: 

PRINT  OUT  OF  GAIN  DATA?  (Y/N)  -> 

REMARKS : 

This  gives  the  user  the  option  to  suppress  the  print 
out  of  the  gain  data  if  it  is  not  needed.  This  saves 
time  and  paper  if  only  the  root  locus  plot  is  desired. 
After  this  step  the  program  runs  to  completion. 


3.  RTLOC  Example 


It  is  desired  to  know  the  root  locus  for  the  open 
loop  transfer  function  G(s)  given  by 


G(s)  = 


1.2  +  s 

s(8  +  9s  +  s2) 


for  gain  K  of  from  0  to  30.  Referring  to  the  steps  des¬ 
cribed  in  section  2,  RTLOC  User's  Guide,  this  problem  would 
be  entered  as  follows: 

STEP  1.  Enter  "EXAMPLE  <CR>"  where  <CR>  =  carriage 

return . 

STEP  2.  Enter  "0,30  <CR>." 

STEP  3.  Enter  "F,l  <CR>." 

STEP  4B.  Enter  "1.2  <CR>  0  <CR>." 

STEP  5.  Enter  "F,3  <CR>  0  <CR>  0  <CR>  1  <CR>  0  <CR> 

8  <CR>  0  <CR> . " 

STEP  6.  Enter  "N  <CR>." 

STEP  9.  Enter  "Y  <CR>." 

The  program  will  then  run  to  completion  producing 
the  output  seen  in  Figure  5. 

The  output  begins  with  the  program  name  and  problem 
identification.  Next  the  transfer  function  numerator  and 
denominator  are  listed  as  coefficients  in  ascending  powers 
of  s  and  then  as  the  real  and  imaginary  parts  of  the  roots 
(i.e.  open  loop  zeros  for  the  numerator  and  open  loop  poles 
for  the  denominator).  The  minimum  and  maximum  values  of  gain 


are  listed  next. 


If  the  option  to  view  only  a  part  of  the  root  locus 
is  selected,  this  fact  is  stated  next  followed  by  the  ranges 
of  sigma  and  omega  that  are  chosen.  (See  Figure  6) 

The  next  several  pages  of  output  contain  the  gain 
data  which  may  be  suppressed  with  the  appropriate  response  in 
step  9.  The  data  point  number  is  listed  with  the  value  of 
gain  and  the  real  and  imaginary  parts  of  the  corresponding 
roots  of  the  open  loop  system. 

The  last  page  of  output  contains  the  root  locus  plot 
itself.  It  begins  with  the  problem  identification  and 
heading  and  ends  with  sufficient  data  listed  to  allow  proper 
interpretation  of  the  plot. 

Figure  6  shows  the  output  from  the  same  example 
problem  with  the  gain  data  suppressed  and  the  option  to 
view  only  a  part  of  the  root  locus  chosen. 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EXAMPLE 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W" 

1.2 

1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

-1.2  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

0 

8 

? 

1 


OPEN-LOOP  POLES 
REAL  PART 


IMAGINARY  PART 


-1 

-8 

0 


0 

0 

0 


MIN  GAIN 
0 


MAX  GAIN 
30 


xxx*xxxxx***xxxxx*xx**xx*****a**xx**x******x**x*** 


Figure  5  RTLOC  Output 
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Figure  5  (Cont.) 
1  GAIN  =  0 


2  GAIN  =  .0492857143 


3  GAIN  =  .  105964286 


4  GAIN  =  . 171144643 


5  GAIN  =  .246102054 


6  GAIN  =  .332303076 


GAIN  =  .431434251 


8  GAIN  =  .545435103 


9  GAIN  =»  .676536033 


10  GAIN  *  .82730221 


ROOTS  ARE 

REAL  PART  I MAG.  PART 

-1  0 

-3  0 

0  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-7.4089 1563E-03  O 
-.998580111  0 

-7.99401098  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.0159695125  0 

-.996917641  0 

-7.98711284  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.0258689696  0 

-.994965313  0 

-7.97916572  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.037328049  0 

-.992664433  0 

-7.97000752  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.0506084588  0 

-.989941249  0 

-7.95945029  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.0660224557  0 

-.986701776  0 

-7.94727577  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.0839456946  0 

-.982824229  0 

-7.93323008  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.104834979  0 

-.978147608  0 

-7.91701742  0 

ROOTS  ARE 

REAL  PART  I MAG.  PART 

-.972453845  0 

-7.89829238  0 

-.129253772  0 


Figure  5 


Figure  5  (Cont.) 

11  GAIN  =  1.60068326 


12  GAIN  *  1.20007146 


13  GAIN  *  1.4293678? 


14  GAIN  =  1.69305879 


15  GAIN  =  1.99630332 


16  GAIN  *  2.34503453 


17  GAIN  =  2.74607543 


18  GAIN  =  3.20727245 


19  GAIN  «  3.73764904 


20  GAIN  *  4.34758211 


ROOTS  ARE 

REAL  PART 
-.965438583 
-7.87665069 
-.  15791073 

ROOTS  ARE 

REAL  PART 
-.956660582 
-7.8516177 
-.  191721725 

ROOTS  ARE 

REAL  PART 
-.945447526 
-7.8226342? 
-.23191818? 

ROOTS  ARE 

REAL  PART 
-7.78903912 
-.28025829? 
-.930702581 

ROOTS  ARE 

REAL  PART 
-7.75004614 
-.339507185 
-.910446682 

ROOTS  ARE 

REAL  PART 
-7.70471555 
-.41482426? 
-.880460137  i 

ROOTS  ARE 

REAL  PART 
-7.65191587 
-.520 147205 
-.827936925 

ROOTS  ARE 

REAL  PART 
-7.59027312 
-.704863441 
-.704863441 

ROOTS  ARE 

REAL  PART 
-7.51310 122 
-.740949392 
-.740949392 

ROOTS  ARE 

REAL  PART 
-7.43330403 
-.78334798? 
-.78334798? 

Figure  5 


I  MAG. 
0 
0 
0 

1MA6 . 
3 
0 
0 

I  MAG. 
0 
0 
0 

IMAG. 

0 

0 

0 

IMAG. 

0 

0 

0 

IMAG. 

0 

0 


PART 


PART 


PART 


PART 


PART 


PART 


IMAG.  PART 


0 

IMAG.  PART 
0 

-.10  1  133243 
.  10  1  133243 

IMAG.  PART 
0 

-.21812379? 

.21312379? 

IMAG.  PART 
0 

-.2970 19384 
.2970 19384 
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Figure 

21 


22 


23 


24 


25 


26 


27 


23 


29 


38 


5  (Cont.) 

GAIN  =  5.84980514 


GAIN  =  5.85564162 


GAIN  =  6.73327358 


GAIN  =  7.35005033 


GAIN  =  9.07684359 


GAIN  =  10.4876559 


GAIN  *  12. 1 100899 


GAIN  =  13.9758891 


GAIN  =  16. 1215582 


GAIN  =  18.5890777 


ROOTS  ARE 

REAL  PART 
-7.33323262 
-.833383692 
-.833333692 

ROOTS  ARE 

REAL  PART 
-7.2144689 
-.392765554 
-.392765554 

ROOTS  ARE 

REAL  PART 
-7.07248055 
-.963759725 
-.963759725 

ROOTS  ARE 

REAL  PART 
-6.90 10354 
-1.0494823 
-1.0494823 

ROOTS  ARE 

REAL  PART 
-6.691 12225 
-1.  15443888 
-1  .  15443888 

ROOTS  ARE 

REAL  PART 
-6.42372813 
-1.28563594 
-1.28563594 

ROOTS  ARE 

REAL  PART 
-6.08944969 
-1.45527515 
-1.45527515 

ROOTS  ARE 

REAL  PART 
-5.62140 192 
-1.68929904 
-1 .63929904 

ROOTS  ARE 

REAL  PART 
-2.6433642 
-4.84656548 
-1.51007033 

ROOTS  ARE 

REAL  PART 
-3.80  149761 
-3.80 149761 
-1.39700479 
Figure  5 


I  MAG.  PART 
0 

-.362332638 

.362882638 

I MAG.  PART 
0 

-.420657402 

.420657402 

I  MAG.  PART 
8 

-.471271626 

.471271626 

I  MAG.  PART 
9 

-.513427903 

.513427903 

I  MAG.  PART 
0 

-.543259781 

.543259731 

I  MAG.  PART 
0 

-.552076334 

.552076334 

I  MAG.  PART 
0 

-.513230319 

.513280319 

I  MAG.  PART 
0 

-.360 139054 
.360 139054 

I  MAG.  PART 
0 
0 
0 

I  MAG.  PART 
-1.23137034 
1.23137034 
0 
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Figure 

31 


32 


33 


5  (Cont .  ) 

GAIN  -  21.426725  ROOTS  ARE 

REAL  PART 
-3.8283818? 
-3.8283818? 
-1.34323622 

GAIN  ■  24 . 6?00 1?5  ROOTS  ARE 

REAL  PART 
-3 . 844?682 
-3 . 844?682 
-1 .3100636 

GAIN  =  28.4428081  ROOTS  ARE 

REAL  PART 
-3.85636161 
-3.85636161 
-1.28727678 


Figure  5 


I  MAG.  PART 
-2. 11737932 
2.  1  1737932 
9 

I  MAG.  PART 
-2 . 7?855?  14 
2 . 7?355?  14 

9 

I  MAG.  PART 
-3.41216555 
3.41216555 
9 


80 


Figure  5  (Cont.) 


PROBLEM  IDENTIFICATION  -  EXAMPLE 

XX  ROOT  LOCUS  PLOT  XX 


ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EXAMPLE  WITH  OPTION 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W" 


1.2 

1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

-1.2  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

0 

3 

? 

1 


OPEN-LOOP  POLES 
REAL  PART 

-1 

-8 

0 


IMAGINARY  PART 

0 

8 

0 


MIN  GAIN  MAX  GAIN 

0  30 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


OPTION  TAKEN 

SIGMA  MIN  *  -2  SIGMA  MAX  =  0 

OMEGA  MIN  *  -1  OMEGA  MAX  =  1 

Figure  6  RTLOC  Output  with  Option 


V.  w' -PLANE  ANALYSIS 


Digital  control  systems  are  becoming  more  and  more  common 
Digital  control  laws  have  unique  characteristics  that  can 
only  be  approximated  by  using  classical  techniques  in  the 
continuous  s  domain. 

In  the  w'  domain  all  analog  control  system  design  tech¬ 
nology  transfers  completely  for  digital  control  system  design 
An  important  advantage  of  the  w'  domain  is  that  non-minimum 
phase  effects  of  the  sampling  and  data-hold  operations  and 
of  sampling  rate  can  be  directly  accounted  for  without 
approximation  while  using  conventional  frequency  domain 
design  and  analysis  tools  such  as  root  locus  and  Bode 
plots.  [Ref.  4] 


A.  BACKGROUND 


This  s  domain  is  used  for  continuous  system  analysis. 
When  a  digital  system  is  considered  the  z,  w,  or  w* domain 
must  be  used.  These  domains  are  related  as  shown  below. 


z-1 

z+1 


sT 

tanh  — 
2 


w'  =  (2/T)  w 


where  T  is  the  sampling  period  and 


tanh  x  = 


x  x 
e  -  e 


x  x 
e  +  e 


In  an  s-plane  root  locus  plot  the  region  of  stability  is 
the  left  half  plane,  that  is  all  roots  with  negative  real 
parts.  This  region  of  stability  is  mapped  into  a  unit  circle 
with  its  center  at  the  origin  in  the  z-plane.  By  use  of 
the  bilateral  transformation  shown  above,  the  stability  region 
of  the  z-plane  is  mapped  back  into  the  left  half  plane  to 
form  the  w-plane.  The  w' -plane  takes  it  a  step  farther  by 
multiplying  the  w-plane  by  the  factor  2/T  where  T  is  the 
sampling  period.  This  gives  the  w' -plane  the  property 
that  not  only  is  the  region  of  stability  the  left  half  plane 
as  in  the  continuous  s-plane  but  w'  approaches  s  as  the  sam¬ 
pling  period  T  approaches  zero  [Ref.  5]. 


TRANSFER  FUNCTION  APPROACH  USING  THE  w' -PLANE 
In  this  section  the  second  order  system 


is  analyzed  using  the  two  transfer  function  programs  dis¬ 
cussed  in  part  IV  of  this  thesis.  It  is  then  converted 
to  the  w' -plane  and  analyzed  for  periods  of  .001  seconds, 

.01  seconds,  and  .1  seconds  using  the  same  two  transfer  func 
tion  programs.  These  results  are  used  to  gain  insight  into 
w' -plane  analysis. 

The  system  is  converted  to  the  w' -plane  by  adding  a 
sampler  and  a  digital  to  analog  converter  in  the  form  of  a 
zero  order  hold.  This  modified  system  is  shown  below. 


-  -sT 
1-e 

s 

K 

_ _ 

i _ s _ 

7 

s( s+10 ) 

The  new  open  loop  transfer  function  is  transformed  to 
the  z-plane  and  then  to  the  w‘ -plane  as  follows: 


G(s)  *  10K( l-e~sT) 


s2(s+10 ) 


G(z)  =  10K(l-z  ) 


‘z 


10 


L<-» 


(l-e-10T)z 

2  10(z-l ) (z-e-10T ) 


=  K 


1-e 


-10T 


z— 1  10(z-e_10T) 


For  the  period 
T=.001 

G(z)  reduces  to 


G(z) 


4.98337  x  10“6  K(z+. 996681) 

-  , 

(z-l)(z-. 990049834) 


Similarly,  for 

T=.01 

4.837418  x  10-4  K(z+. 9672185) 

G(z)  =  - 

(z-l)(z-. 904837418) 


And  for 


T= .  1 

.03678794412  K( z+ . 718281827 ) 

G(z)  =  - 

(z-l)(z-. 3678794412) 

Now  G(z)  is  converted  to  the  w'  domain  using  the  rela 


tionship 


1  +  T/2  w' 

z  =  - 

1  -  T/2  w' 

For  period 

T=. 001 

1  +  .0005  w' 

z  =  - 

1  -  .0005  w* 


Substituting  into  equation  (1)  gives 


4.98337  x  10“6  K 


G(w')  = 


1+.0005W* 

-  +  .996681 

l-.0005w'  J 


l+.0005w' 

1-.0005W* 


-  1 


1+.0005W* 

l-.0005w' 


-  .990049834 


which  reduces  to 


-4.156878  x  10“9  K  (w*  +  1202819)  (w'  -  2000) 

G(w')  -  -  ( 

w'  (w*  +  9.9999) 

Note  that  the  gain  is  1.0K  which  is  the  same  as  the  s  domain 
gain.  Similarly,  for 

T=.01 


1+.005W' 

z=  - 

l-.005w’ 
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-4.1625027  x  10  K(w • +12002 . 00418 ) (w 1 -200 ) 

G(wl)  =  -  (3) 

w'  (w*  +  9.991674985) 

And  for 

T=.  1 


l+.05w' 

z  -  - 

l-.05w' 

-.0037882843  K(w*  +  121 . 9858708 ) (w‘ -20 ) 

G(W»)  =  -  (4) 

w'  (w*  +  9.242343139) 

Note  again  that  the  gain  is  1.0K  also  in  equations  (3)  and 
(4). 

1.  Frequency  Response 
a.  s-Plane 

The  open  loop  transfer  function  of  the  system 
in  the  s  domain  is 

10K 

G(s)  =  - 

s(s+10 ) 

The  frequency  response  of  G(s)  with  a  gain  K  of  5  is  described 
by  the  output  of  the  FRESP  program  shown  in  Figure  7.  From 
the  Bode  plots  it  is  seen  that  the  gain  margin  is  infinite 
and  the  phase  margin  is  65.1  degrees. 


b.  w' -Plane  with  a  Period  of  .001  Seconds 

The  open  loop  transfer  function  for  the  equiva¬ 
lent  sampled  system  is  represented  by  equation  (2)  for  a 
period  T  of  .001  seconds.  This  transfer  function  is 
entered  into  the  FRESP  program  as  if  w'  were  an  s.  The 
resulting  output  in  Figure  8  is  interpreted  as  described 
earlier  in  the  thesis  keeping  in  mind  that  all  references 
to  s  are  actually  referring  to  w' .  From  this  output  it 
is  seen  that  the  gain  margin  is  no  longer  infinite  as  it 
was  in  the  continuous  case  but  it  is  still  quite  high  at  50 
dB.  The  phase  margin  has  dropped  slightly  from  65.1  degrees 
to  64.6  degrees. 

c.  w' -Plane  with  a  Period  of  .01  Seconds 

If  the  sampling  rate  of  the  system  is  decreased 
so  that  the  sampling  period  is  .01  seconds,  then  the  open 
loop  transfer  function  is  represented  by  equation  (3). 
Inputting  the  transfer  function  into  FRESP  results  in  the 
output  shown  in  Figure  9.  From  the  Bode  plots  the  gain 
margin  is  found  to  be  32  dB.  This  is  down  from  50  dB  for 
the  case  of  a  period  of  .001  seconds  indicating  a  decrease 
in  stability.  The  phase  margin  is  64.5  degrees  which  is  also 
lower  than  the  previous  case  but  only  by  one  tenth  of  a 
degree. 

d.  w'-Plane  with  a  Period  of  .1  Seconds 

If  the  period  is  again  increased  to  .1  seconds, 
the  open  loop  transfer  function  is  represented  by  equation 


(4).  Entering  this  transfer  function  into  FRESP  produces 
the  output  seen  in  Figure  10.  From  this  output  the  gain 
margin  is  found  to  have  decreased  to  13.2  dB  and  the  phase 
margin  to  53.7  degrees. 

e.  Summary  of  Frequency  Response  Results 

The  table  below  is  a  brief  summary  of  the  gain 
and  phase  margins  found  in  each  case. 


CASE 


GAIN  MARGIN 


PHASE  MARGIN 


s-plane 
w'  ,  T=.001 
w'  ,  T= . 01 
w ' ,  T=.l 


infinite 
50  dB 
32  dB 
13.2  dB 


65.1  deg. 

64.6  deg. 
64.5  deg. 

53.7  deg. 


It  can  be  seen  from  this  table  that  the  con¬ 
tinuous  case  is  the  most  stable.  The  sampled  cases  become 
less  stable  as  the  sampling  period  increases.  It  is  also 
noticed  that  the  gain  margin  is  more  sensitive  to  changes 
in  the  sampling  period  than  the  phase  margin. 

2 .  Root  Locus 

Now  a  comparison  is  made  using  the  root  locus  program 
using  the  same  cases  used  above, 
a.  s-Plane 

If  the  open  loop  transfer  function 


10K 

G(s)  =  - 

s(s+10) 


is  entered  into  the  RTLOC  program,  the  result  is  the  output 
in  Figure  11.  The  root  locus  plot  shows  that  the  system 
never  becomes  unstable  at  any  gain. 
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b.  w' -Plane  with  a  Period  of  .001  Seconds 

If  the  equivalent  sampled  transfer  function  with  a 
sampling  period  of  .001  seconds  represented  by  equation  (2) 
is  entered  into  RTLOC,  the  output  will  be  that  found  in 
Figure  12.  From  the  root  locus  plot  it  can  be  seen  that 
there  is  a  slight  tendency  for  the  plot  to  curve  toward  the 
right  half  plane  as  it  moves  further  from  the  real  axis.  T 
tendency  is  so  slight  that  it  will  take  a  very  large  gain 
to  drive  the  system  unstable. 

c.  w' -Plane  with  a  Period  of  .01  Seconds 

If  the  sampling  period  is  increased  to  .01  seconds 
the  transfer  function  is  represented  by  equation  (3). 
Entering  this  transfer  function  into  RTLOC  results  in  the 
output  seen  in  Figure  13.  It  can  be  seen  that  the  tendency 
for  the  curve  to  bend  toward  the  unstable  right  half  plane 
is  increased  indicating  that  a  lower  value  of  gain  than  in 
the  previous  case  will  drive  the  system  to  instability. 

d.  w' -Plane  with  a  Period  of  .1  Seconds 

When  the  sampling  period  is  further  increased  to  .1 
seconds  represented  by  equation  (4),  the  resulting  RTLOC 
output  is  that  seen  in  Figure  14.  In  the  root  locus  plot 
for  this  case  the  tendency  to  become  unstable  is  much  more 
pronounced.  Here  the  plot  bends  into  the  right  half  plane 
within  the  limited  boundaries  of  the  portion  plotted.  It 
becomes  unstable  for  values  of  gain  K  greater  then  24.6. 


e.  Summary  of  Root  Locus  Comparison 

It  can  be  seen  that  the  continuous  system 

10K 

G(s)  =  - 

s(s+10) 

is  always  stable.  When  this  system  is  sampled  it  can  be 
seen  by  equations  (2),  (3),  and  (4)  that  a  zero  is  added 
in  the  right  half  plane  making  it  a  non-minimum  phase 
system.  It  can  also  be  seen  that  by  decreasing  or  increasing 
the  sampling  period  the  distance  of  this  zero  from  the 
origin  increases  or  decreases  respectively.  The  closer 
this  zero  is  to  the  imaginary  axis,  the  greater  effect  it 
has  on  bending  the  root  locus  into  the  right  half  plane. 

This  effect  can  be  seen  by  examining  the  root  locus  plots 
in  Figures  11  through  14. 
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Figure  8  w' -Plane  Frequency  Response  Example,  T=.001 
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-2.84914171 

2.1888924 

-1.889367 

-183.66913 

2.45488781 

-.498358539 

-1.75212349 

1.81944588 

-1.84367592 

-185.434887 

3.83385554 

-.488148485 

-1.49852834 

1.54595572 

-1.88243386 

-187.855588 

3.51119177 

-.447294883 

-1.25981191 

1.34348584 

-1.92598991 

-118.351133 

4.8378173 

-.451325898 

-1.85419372 

1.14858184 

-1.97442499 

-113.137719 

4.44138888 

-.431817282 

-.874428211 

.97721191 

-2.82849743 

-116.224383 

3.33449929 

-.488474439 

-.718744596 

.826718149 

-2.88759424 

-119.418383 

4.13598734 

-.381235751 

-.388778541 

.494719472 

-2.15147917 

-123.28218 

7.83488239 

-.338349774 

-.441388989 

.579324934 

-2.22825484 

-127.211277 

8.1113884 

-.314458228 

-.359537358 

.47897869 

-2.29255364 

-131.353694 

9.32483358 

-.288578444 

-.274244392 

.392368887 

-2.34754996 

-135.651815 

18.7224724 

-.244887234 

-.28444584 

.318348855 

-2.444138 

-148.838842 

12.3284474 

-.288143943 

-.148747473 

.255843176 

-2.5218441 

-144.445353 

14.1747418 

-.174282984 

-.185522173 

.283738772 

-2.59714482 

-148.886636 

14.2975884 

-.143437538 

-.8728478818 

.168884843 

-2.47154189 

-153.86923 

18.7381745 

-.114242424 

-.8488738445 

.126899635 

-2.74357414 

-157.19539 

21.5443472 

-.8929483537 

-.8317858744 

.8982874138 

-2.81285749 

-161.164932 

24.7787439 

-.8734841787 

-.8197248894 

.8748854428 

-2.87935849 

-164.974781 

28.4883591 

-.8573538451 

-.8115449167 

.8587844732 

-2.94325927 

-168.636395 

32.7434921 

-.8447358414 

-4.14914418E-83 

.8451564981 

-3.88499418 

-172.173544 

37.4493587 

-.8345424152 

-2.64835976E-83 

.834663932 

-3.84511716 

-175.41834 

43.2874135 

-.8245739352 

-4.S9917489E-G4 

.8265799144 

-3.12428864 

-179.888617 

49.7782345 

-.8283593343 

8.48884793E-84 

.8283769894 

3.89994478 

177.414942 

37.2234774 

-.8135327833 

1.57682436E-83 

.8156324331 

3.84855195 

174.218856 

45.7933234 

-.8118553889 

1.9329344E-83 

.8128119388 

2.97997212 

178.739887 

73.4443342 

-9.82245397E-83 

2.83532664E-83 

9.2537932E-83 

2.9176184 

147.147281 

84.9749819 

•4.8589325 IE-83 

2.834 1848E-83 

7.15478944E-83 

2.85381322 

163.465675 

188.888882 

-5.21884889E-83 

1. 93572 133E-83 

S.55883487E-83 

2.78584384 

159.418253 

Figure  9 
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ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 

MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =*  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  =  100  RADIANS/SEC 
AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-40  DECIBELS 


Figure  9 
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ABSCISSA  ->  corn  ON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  90  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  100  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  -  180  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  -  -188  DEGREES 


Fig1  re  9 
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FREQUENCY  RESPONSE 

PROBLEM  IDENTIFICATION  -  EX  1  W' -PLANE  T= . 1 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXX 


GAIN*-. 0189414215 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 


-2439.71742 

101.985871 

1 


NUMERATOR  ROOTS  ARE 
REAL  PART 


IMAGINARY  PART 


-121.985871  9 

20  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 
0 

9.24234314 
1  ' 

DENOMINATOR  ROOTS  ARE 
REAL  PART 


0 


IMAGINARY  PART 

0 

0 


-9.24234314 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


Figure  10 


w‘ -Plane  Frequency  Response  Example,  T=.l 
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Figure  10  (Cont.  VrqbLEM  IDENTIFICATION  -  EX1  W'-PLANE  T= .  1 


RADIAN  FREQ. 

REAL  PART 

. 117224722 

-.749981417 

.123284474 

-.749849949 

.141747414 

-.749828888 

.142973184 

-.749772739 

.187381743 

-.749499429 

.213443449 

-.749482979 

.247717434 

-.749473252 

.284813388 

-.749384472 

.327454917 

-.749883474 

.374493582 

-.748788894 

.43287413 

-.748399844 

.497712359 

-.747884175 

.572234749 

-.747288244 

.457933228 

-.744313994 

.734443332 

-.745135214 

.849749808 

-.743582787 

1.08888801 

-.741548723 

1.149757 

-.738858839 

1.32194114 

-.735343947 

1.51991109 

-.730749803 

1.74752841 

-.724744284 

2.88923382 

-.717809119 

2.31812972 

-.787012534 

2.45488781 

-.494228885 

3.83385554 

-.478841194 

3.51119177 

-.457794328 

4.8378173 

-.432848475 

4.44158888 

-.482734245 

5.33449929 

-.547282992 

4.13598734 

-.524334324 

7.05488239 

-.488842854 

8.1113884 

-.431914354 

9.32483358 

-.38117579 

10.7224724 

-.338541899 

12.3284474 

-.281974491 

14.1747418 

-.237827582 

14.2975884 

-.194848811 

18.7381745 

-.141992858 

21.3443472 

-.132349482 

24.7707439 

-.10827854 

28.4883591 

-.0883944878 

32.7434921 

-.8728831383 

37.4493587 

-.848493815 

43.2874133 

-.8388148715 

49.7782343 

-.8433112738 

57.2234774 

-.8373272873 

45.7933234 

-.8338884779 

73.4443342 

-.8294938247 

84.9749819 

-.8271844397 

IN  .888882 

-.0231332985 

IMAGINARY  PART 

MAGNITUDE 

-44.4214872 

44.4277179 

-40.5447918 

48.5537253 

-35.2428823 

35.2707734 

-38.4444494 

38.4758137 

-24.4484887 

24.4792143 

-23.1989122 

23.2838238 

-28.145587 

28.1794298 

-17.5334525 

17.5494543 

-15.2434888 

15.2418832 

-13.2507881 

13.271848 

-11.5144835 

11.5487753 

-18.9849112 

18.8348194 

-8.4925514 

8.72448734 

-7.54777487 

7.58458359 

-4.55824855 

4.59251428 

-5.48059274 

5.72985309 

-4.92181445 

4.97734475 

-4.25918454 

4.32279798 

-3.4798489 

3.75248159 

-3.17248843 

3.25547824 

-2.72772751 

2.82237198 

-2.33475582 

2.44428514 

-1.99239481 

2.11412199 

-1.48848289 

1.82555418 

-1.41949479 

1.57311942 

-1.1813151 

1.35211007 

-.978377373 

1.15851483 

-.784821224 

.988939817 

-.42833237 

.848554272 

-.478807285 

.71181473 

-.354144355 

.598388931 

-.253984373 

.581059494 

-.178419147 

.417419384 

-.18474951 

.344747474 

-.8544999148 

.287233844 

-.8182471584 

.237738343 

4.90255474E-03 

.194948999 

.9231913419 

.14344458? 

.8327938441 

.1345451? 

.8375858315 

.114414214 

.8398773822 

.8948299544 

.9384219724 

.8823985323 

.9344598583 

.878438582 

.9337744111 

.0418144784 

.8387785931 

.0531298584 

.8274948923 

.0444412814 

.8247128748 

.8412985915 

.921988753 

.0349814454 

.8193382734 

.8332913884 

.8149944478 

.8383397445 

Figure  10 
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PHASE(RAO) 

PHASE(OEG) 

-1.58487981 

-98.9215483 

-1.58928821 

-91.0595392 

-1.59285717 

-91.213189 

-1.59524843 

-91.4085381 

-1.59898059 

-91.4102882 

-1.48318824 

-91.8513782 

-1.48794548 

-92.1285223 

-1.41358421 

-92.4471289 

-1.41989834 

-92.8133714 

-1.42724549 

-93.2343435 

-1.43549823 

-93.7181881 

-1.44539445 

-94.274203 

-1.45454519 

-94.9138817 

-1.44935447 

-95.4478111 

-1.4848459 

-94.4899828 

-1.78895514 

-97.4575854 

-1.72833548 

-98.5488892 

-1.74254117 

-99.8414341 

-1.7488284 

-181.388402 

-1.79717901 

-182.978889 

-1.83849917 

-184.379914 

-1.84851412 

-187.858124 

-1.91179885 

-109.53754 

-1.94089918 

-112.351287 

-2.01441584 

-115.532159 

-2.87887135 

-119.110597 

-2.14878819 

-123.111955 

-2.22422829 

-127.553873 

-2.31148514 

-132.438391 

-2.48429439 

-137.754085 

-2.58418884 

-143.474919 

-2.41888982 

-149.5424 

-2.72872574 

-155.884158 

-2.83444442 

-142.414377 

-2.94998582 

-149.821752 

-3.8444771 

-175.593124 

3.18454824 

177.991788 

2.99939485 

171.852794 

2.89989541 

144.184802 

2.887498 94 

140.85744 

2.72418143 

154.198744 

2.45444429 

152.2831 

2.59919124 

148.922743 

2.55495334 

144.388895 

2.5238842 

144.488879 

2.58578445 

143.578934 

2.50809415 

143.245809 

2.58593127 

143.579337 

2.52287784 

144.584447 

2.547031? 

145.93423 

ABSCISSA  ->  COMMON  LOO  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 

MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  »  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOW  ON  ABSCISSA  =*  100  RADIANS/SEC 
AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-40  DECIBELS 


Figure  10 


ABSCISSA  ->  COt“W ON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  90  DEGREES 
MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOW  ON  ABSCISSA  -  100  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  -  180  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  -  -180  DEGREES 


Figure  10 


ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EX  1  S-PLANE 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W" 


1 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 

0 

10 

1 


OPEN-LOOP  POLES 
REAL  PART 

0 

-10 

MIN  GAIN 
0 


IMAGINARY  PART 

0 

0 

MAX  GAIN 
300 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


Figure  11  s-Plane  Root  Locus  Example 
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Figure  11  (Cont.) 


1 

GAIN 

3 

0 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-10 

0 

0 

0 

2 

GAIN 

» 

.492857143 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.0495310467 

0 

-9.95046896 

0 

3 

GAIN 

SS 

1.05964286 

ROOTS 

ARE 

REAL  PART 

IMA6. 

PART 

-. 107111575 

0 

-9.89288843 

0 

4 

GAIN 

= 

1.71144643 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.  174178456 

0 

-9.82582155 

0 

5 

GAIN 

3 

2.46102054 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.252476491 

0 

-9.74752351 

0 

6 

GAIN 

3 

3.32303076 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.344146776 

0 

-9.65585323 

0 

7 

GAIN 

= 

4.31434252 

ROOTS 

ARE 

REAL  PART 

I  MAG . 

PART 

-.451851202 

0 

-9.5481488 

0 

8 

GAIN 

3 

5.45435104 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.578953861 

0 

• 

-9.42104614 

0 

9 

GAIN 

3 

6.76536084 

ROOTS 

ARE 

REAL  PART 

I  MAG. 

PART 

-.729796355 

0 

-9.27028365 

0 

10 

GAIN 

* 

8.27302211 

ROOTS 

ARE 

REAL  PART 

IMAG. 

PART 

-.910 137178 

0 

-9.08986283 

0 

1  1 

GAIN 

m 

10.0068326 

ROOTS 

ARE 

REAL  PART 

IMAG. 

PART 

-1. 12789883 

0 

-8.87210117 

0 

12 

GAIN 

3 

12.0007146 

ROOTS 

ARE 

REAL  PART 

IMAG. 

PART 

-1.39454782 

0 

-8.60545218 

0 

13 

GAIN 

* 

14.2936789 

ROOTS 

ARE 

REAL  PART 

IMAG. 

PART 

-1.72794849 

0 

-8.27205151 

0 

Figure  11 
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Figure  11  (Cont.) 


14 

GAIN 

a 

16.9305879 

ROOTS. 

ARE 

REAL 

PART 

I MAG.  PART 

-2. 15932893 

0 

-7.84067107 

0 

15 

GAIN 

a 

19.9630332 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-2.75568122 

0 

-7.24431878 

0 

16 

GAIN 

a 

23.4503453 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-3.75514874 

0 

-6.24485127 

0 

17 

GAIN 

a 

27.4607543 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-1 .56867915 

-5 

1 .56867915 

18 

GAIN 

a 

32.0727246 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-2.65945945 

-5 

2.65945945 

19 

GAIN 

a 

37.3764904 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-3.51802365 

-5 

3.51802365 

28 

GAIN 

a 

43.4758212 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-4.29835098 

-5 

4.29835098 

21 

GAIN 

a 

50.4900515 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-5.04876732 

-5 

5.04876732 

22 

GAIN 

a 

58.5564163 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-5.79279003 

-5 

5.79279003 

23 

GAIN 

a 

67.8327359 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-6.54467233 

-5 

6.54467233 

24 

GAIN 

a 

78.5005034 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-7.31440383 

-5 

7.31440383 

25 

GAIN 

a 

90.7684361 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-8. 10977411 

-5 

8.  10977411 

26 

GAIN 

a 

104.876559 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-8.93736867 

-5 

8.93736867 

Figure 
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Figure  11  (Cont.) 

27  GAIN  -  121. 1009 


28  GAIN  =  139.758892 


29  GAIN  »  161.215582 


30  GAIN  *  185.890777 


31  GAIN  »  214.267251 


32  GAIN  -  246.900195 


33  GAIN  -  284.428082 


ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-9.30318663 

-5 

9.80310663 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-10 .7125577 

-5 

10.7125577 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-11.6711431 

-5 

11.6711431 

ROOTS 

ARE 

REAL 

PART 

I MAG.  PART 

-5 

-12.6842728 

-5 

12.6842728 
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REAL 

PART 

I MAG.  PART 
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-13.7574435 
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13.7574435 
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-5 
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-5 

16. 1067713 
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Figure  11  (Cont.) 


v. 


PROBLEM  IDENTIFICATION  -  EX  I  S-PLANE 
XX  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  (SIGMA)  AXIS 
OROINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  l 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -22  TO  II 
ORDINATE,  -17  TO  16 


ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  EX  1  W'-PLANE  T=.001 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W' 

-2.405638E+09 

1200819 

1 


OPEN-LOOP  ZEROS 

REAL  PART  IMAGINARY  PART 

-1202819  0 

2000  0 


DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W 
0 

9.9999 

1 


OPEN-LOOP  POLES 
REAL  PART 

0 

-9.9999 


IMAGINARY  PART 

0 

0 


MIN  GAIN 
0 


MAX  GAIN 
-1.2346E-07 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


Figure  12  w' -Plane  Root  Locus  Example,  T=.001 


Figure  12  (Cont.) 
1  GAIN  =  8 
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REAL  PART 
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-9.9999 
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ROOTS 

ARE 

REAL  PART 

I MAG.  PART 

-.346367788 

0 

-9.65186347 

0 

GAIN 

a 

-1 .80425804E-09 
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I MAG.  PART 
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-.734883239 

0 

-9.26161935 

0 

GAIN 

a 

-3 . 45977784E-09 
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I MAG.  PART 

-.916725984 
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I MAG.  PART 
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I MAG.  PART 
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Figure  13  w' -Plane  Root  Locus  Example,  T=.01 
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Figure  13  (Cont.) 
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Figure  14  (Coat.) 
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Figure  14  (Cont.) 
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C.  ROOT  LOCUS  TEMPLATES  IN  w* -PLANE 

The  s-plane  consists  of  an  abscissa  representing  the 
real  portion  of  the  complex  variable  s  and  an  ordinate 
representing  the  imaginary  part  where 
s=  G  +  3  Ci) 

The  w* -plane  approaches  the  s-plane  as  the  sampling 
period  T  approaches  zero.  (See  [Ref.  11]  for  a  proof.) 

The  w' -plane,  like  the  s-plane,  consists  of  an  abscissa 
representing  the  real  part  and  an  ordinate  representing 
the  imaginary  part  where 
w'=u'+j  V' 

The  s  and  w'  planes  are  related  as  shown  below. 

u'  =  2/T  tanh  (  (J  T/2) 

V'  ■  2/T  tan(  (J  • T/2) 

O’  =  2/T  tanh-1  (u'T/2)  where  -1  <  u'T/2  <  1 

GJj  =  2/5  tan-1  (  ^ ' T/2)  where  -1  <  ^'T/2  <  1 

In  the  s-plane,  O’  represents  damping.  Therefore  a  line 
of  constant  damping  in  the  s-plane  is  a  vertical  line 
parallel  to  the  imaginary  axis  and  passing  through  the  proper 
value  of  O’  on  the  real  axis.  Since  the  real  axis  in  the 
w' -plane  represents  u'  which  is  not  the  damping  but  only 

related  to  O’  as  shown  above,  it  is  not  possible  to  easily 

determine  the  damping  of  a  particular  root  by  simply  ob¬ 
serving  the  position  of  the  root  in  the  plane.  This  is  also 
true  of  the  damped  natural  frequency,  the  natural  frequency, 
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and  the  damping  ratio.  Because  of  this  it  is  helpful  to 
have  templates  of  constant  parameters  in  the  w' -plane  to 
make  interpretation  of  the  characteristics  of  a  root 
by  its  location  easier.  The  following  sections  give  some 
insight  into  the  nature  of  these  templates  in  the  w' -plane. 

1.  Constant  Damping 

The  constant  damping  templates  are  shown  in  Figure  15 
There  are  three  templates,  one  of  each  of  three  values  of 
sampling  period,  .001  seconds,  .01  seconds,  and  .1  seconds. 
Each  template  shows  lines  of  constant  damping  for  the  values 
25,  50,  75,  and  100.  These  templates  are  created  by  using 
the  above  relationship  between  (J  and  u'  to  find  the  constant 
value  of  u'  that  corresponds  to  the  constant  value  of  (J  and 
then  plotting  the  constant  u’  value.  From  the  templates  it 
is  seen  that  for  a  period  of  .001  there  is  practically  no 
difference  from  the  s-plane.  For  a  period  of  .01  seconds, 
some  distortion  is  noticed  as  the  values  of  damping  approach 
the  reciprocal  of  the  value  of  the  period.  This  distortion 
is  in  the  form  of  u'  becoming  smaller  than  the  value  of 
damping.  And  for  a  period  of  .1  seconds  gross  distortion 
is  seen.  Figure  15  shows  an  extra  template  for  a  period 
of  .1  seconds  but  for  values  of  damping  of  1,  2,  3,  and  4. 

In  this  case  only  little  distortion  is  present  because  the 
values  of  damping  are  small  compared  to  the  reciprocal  of 


the  period. 


2 .  Constant  Damped  Natural  Frequency 
Figure  16  shows  the  templates  for  constant  damped 

natural  frequency  for  values  of  damped  natural  frequency  of 
25,  50,  75,  and  100.  The  results  here  are  similar  to  the 
constant  damping  case.  As  the  period  increases,  the 
distortion  between  the  value  of  Cdi  and  the  value  of  in¬ 
creases.  In  this  case  the  value  of  V1  becomes  greater  than 
the  value  of  damped  natural  frequency  as  it  distorts. 

In  the  case  where  the  period  is  .1  seconds  the  distortion 
is  so  great  that  the  template  is  useless  and  negative 
values  of  J/1  are  produced  by  the  tangent  function  in  the 
relationship. 

3.  Constant  Natural  Frequency 
Templates  for  constant  natural  frequency  in  the 

w' -plane  are  shown  in  Figure  17.  In  the  s-plane  constant 
natural  frequency  C)Jn  plots  as  a  circle  with  its  center  at 
the  origin  where 

=  sqrt  (  CT  2  +  CJj 2  ) 

To  plot  constant  natural  frequency  in  the  w' -plane,  substitute 

u'  =  2/T  tanh(  <J  T/2) 

V'  -  2/T  tan(  U*  T/2) 

into  the  above  equation  and  solve  for  V'  to  get 

V'  =  2/T  tan  [  T/2  (CJ2-  [2/T  tanh-1  (u’T/2)]2}* 
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For  a  constant  value  of  natural  frequency  and  period  this 


equation  is  used  to  plot  u'  vs  .  In  the  w'-plane  for  a 
period  of  .001  seconds  and  a  value  of  natural  frequency  of  25 
it  plots  very  close  to  a  perfect  circle.  As  the  values  of 
natural  frequency  increase  they  form  concentric  ellipses 
which  are  elongated  along  the  imaginary  axis  and  contracted 
along  the  real  axis.  This  effect  becomes  more  dramatic 
as  the  period  increases.  An  extra  template  for  a  period  of 
.05  seconds  is  added  to  better  show  the  trend  since  the 
change  from  .01  seconds  to  .1  seconds  is  so  great. 

4.  Constant  Damping  Ratio 

The  damping  ratio,  £  ,  can  be  defined  as 

£  =  |  (T  /  CJn  |  ■  cos  Q 

where  J  |  represents  the  absolute  value  and  Q  is  the  angle 
formed  by  the  negative  real  axis  and  a  line  joining  the 
origin  and  the  root  in  a  root  locus  plot.  In  the  s-plane  a 
plot  of  constant  damping  ratio  is  a  line  radially  out  from 
the  origin  forming  the  correct  angle  $  with  the  negative 
real  axis.  To  plot  constant  damping  ratios  in  the  w’-plane, 
substitute  the  proper  relations  into  the  above  equation  to 
get 

^  tanh-1  (u’T/2) 


(u'T/2)]  +  [tan-x  {  V'  T/2)]^} 


Solving  for  u'  gives 


1 

tan"1  (  V1  T/2 ) £ 

u'  =  2/T  tanh 

sqrt ( 1-  ) 

This  equation  is  used  to  plot  constant  damping  ratio  for 
a  given  sampling  period.  Figure  18  shows  templates  for 
constant  damping  ratio  for  values  of  damping  ratio  of  .1, 

.5,  .707,  and  .9.  From  these  templates  it  can  be  seen  that 
there  is  a  region  near  the  origin  that  can  be  interpreted 
the  same  as  the  s-plane  for  a  close  approximation  and  that 
the  size  of  this  region- depends  on  the  value  of  the  sampling 
period.  Beyond  this  region  a  template  makes  the  interpreta¬ 
tion  easier. 
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CONSTANT  DAMPING  W' -PLANE  T=.001 

X*  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  <U>  AXIS 

Ul  ORDINATE  ->  IMAGINARY  <NU>  AXIS 

M  TIC  MARKS  SHOW  INTERVALS  OF  5 

■  THE  PLOT  FRAME  LIMITS  ARE! 

ABSCISSA,  -163  TO  37 
ORDINATE,  -180  TO  100 


Figure  15  Constant  Damping  Templates 
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ABSCISSA  ->  REAL  <U>  AXIS 
ORDINATE  ->  IMAGINARY  <NU>  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  S 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -160  TO  40 
ORDINATE,  -100  TO  100 


Figure  15 


ABSCISSA  ->  REAL  <U>  AXIS 
ORDINATE  ->  IMAGINARY  <NU>  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  5 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -118  TO  82 
ORDINATE,  -100  TO  100 


Figure  15 


Figure  15 


CONSTANT  DAMPING  W' -PLANE  T= . 1 

XX  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  <U>  AXIS 
ORDINATE  ->  IMAGINARY  <NU>  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  1 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -7  TO  3 
ORDINATE,  -5  TO  5 


Figure  15 


ABSCISSA  ->  REAL  < U)  AXIS 
ORDINATE  ->  IMAGINARY  (NU)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  5 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -190  TO  100 
ORDINATE,  -37  TO  163 


Figure  16  Constant  Damped  Natural  Frequency  Templates 
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Figure  16  (Coat.) 
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CONSTANT  NATURAL  FREQUENCY  W' -PLANE  T-.081 
XX  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  <U)  AXIS 
ORDINATE  ->  IMAGINARY  <NU>  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  5 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -128  TO  77 
ORDINATE,  - 185  TO  100 


fjP  Figure  17 


Figure  17  (Cont.) 


CONSTANT  NATURAL  FREQUENCY  W' -PLANE  T-.i 
XX  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  <U>  AXIS 
ORDINATE  ->  IMAGINARY  CNU)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  S 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -120  TO  85 
ORDINATE,  -105  TO  100 


Figure  17 
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CONSTANT  DAMPING  RATIO  W' -PLANE  T=.80t 
XX  ROOT  LOCUS  PLOT  XX 


THE  PLOT  FRAME  LIMITS  ARE: 

ABSCISSA,  -3800  TO  1008 
ORDINATE,  0  TO  4000 
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Figure  18  Constant  Damping  Ratio  Templates 
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18  (Cont . ) 


CONSTANT  DAMPING  RATIO  W' -PLANE  T=.0 1 
XX  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  <U>  AXIS 
ORDINATE  ->  IMAGINARY  <NU>  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  10 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -300  TO  100 
ORDINATE,  0  TO  400 


Figure  18 
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Figure  18  (Cont.) 
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CONSTANT  DAMPING  RATIO  W' -PLANE  T-. 1 
XX  ROOT  LOCUS  PLOT  XX 


ABSCISSA  ->  REAL  <U>  AXIS 
ORDINATE  ->  IMAGINARY  <NU>  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  l 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -30  TO  10 
ORDINATE,  0  TO  40 


Figure  18 


Converting  to  the  z  domain  gives 


9  (z+1) 

G-  (z)  =  K9T^  - 9 

1  2  (z-1)2 


where 


K2  = 


KjKt 

2J 


59.8 


Converting  this  to  the  w  domain  give 


KqT2  1-w 
G1<W)  =  -A—  _ 

2  w 


where 


w  =  w'T/2 


The  open  loop  transfer  function  is 
G(w)  *  KpK^Cw) 


where 


this  transfer  function  W!ff§  entered  into  the  FRESP  program. 

The  resulting  Bode  plots  in  Figure  20  show  that  the  compen¬ 
sated  system  is  stable  with  a  gain  margin  of  10.5  dB  and 
phase  margin  of  41  degrees. 

Figure  21  shows  the  output  from  the  RTLOC  program  for  the 
compensated  system.  This  figure  includes  the  normal  root 
locus  in  addition  to  an  expanded  portion  of  the  root  locus 
to  obtain  more  detail.  These  plots  are  of  w  and  not  w'  in 
this  case.  Values  of  w'  can  be  found  by  the  simple  relation 

w'  =  2w/T 
where 

T  =  .01  seconds 
for  this  problem. 

On  the  expanded  root  locus  plot  in  Figure  21,  a  root 
associated  with  a  certain  gain  may  be  identified  and 
templates  similar  to  the  ones  discussed  in  section  V.C.  may 
be  used  to  easily  obtain  the  characteristics  associated  with 


the  root 
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Figure  19  Uncompensated  System  Frequency  Response 
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PROBLEM  I  DENT I F I  CAT I  ON 

-  KATZ 

D<  AMPLE 

RADIAN  FREQ. 

REAL  PART 

IMAGINARY  PART 

MAGNITUDE 

PHASE(RAD) 

PHASE (DEG) 

m 

r 

.184761575 

-2.88455884 

.21099866 

2.81552885 

3.83721188 

174.319434 

\ 

w 

.1149757 

-1.44421932 

.19134478 

1.4751832 

3.82711946 

173.441242 

.126185688 

-1.38166411 

.174346238 

1.39262867 

3.91607043 

172.308168 

K*i% 

.138488637 

-1.14788182 

.158857798 

1.15882958 

3.00397938 

172.115482 

• 

-y* 

.151991119 

-.952327477 

.144745389 

.963264672 

2.99875683 

171.357759 

is 

.166818854 

-.798639881 

.131886535 

.381563527 

2.97638452 

178.529748 

■  i  ms 

.183873829 

-.656492388 

.128178898 

.667311733 

2.968524 

169.425591 

m 

.288923381 

-.544954794 

.189494518 

.555347942 

2.94330951 

168.639273 

r's 

.228513875 

-.452432784 

.8997473267 

.463382138 

2.92455384 

167.544407 

»";n 

.242812828 

-.375617579 

.8999842724 

.386461962 

2.90414539 

166.295333 

►V* 

.26568878 

-.311844312 

.9828285873 

.322656861 

2.38197822 

165.125248 

.291585388 

-.258898627 

.0754783239 

.249474375 

2.35794731 

163.743377 

'T 

.319926714 

-.214942188 

.9687657482 

.225674241 

2.83195624 

162.259198 

ias 

.351119174 

-.17844878 

.9626547886 

.139129143 

2.38392119 

168.652908 

**> 

.385352862 

-.148151312 

.8578985323 

.158778714 

2.77377459 

158.925748 

•/ 

.422924291 

-.122997822 

.0528187467 

.133545559 

2.74148151 

157.075374 

.V 

.464158887 

-.182114953 

.0473975628 

.112578328 

2.78782483 

155.181268 

.589413886 

-.9847774274 

.9431868939 

.8951438593 

2.47044244 

153.995136 

* 

.559881823 

-.8783838757 

.0393502893 

.9886370586 

2.63188424 

158.791329 

.613598733 

-.8584339381 

.035854518 

.0485570613 

2.5912399 

143.467143 

*■* « 

.473415872 

-.8485128753 

.9326693915 

.8584874545 

2.54893266 

146.943136 

.739872211 

-.8482742413 

.8297678587 

.8508824612 

2.58512214 

143.532977 

3 

.31113884 

-.8334388431 

.027122628 

.8438558773 

2.46818142 

148.953479 

» 

.898215895 

-.8277488588 

.8247131284 

.0371472377 

2.41421 

133.324893 

.977889969 

-.923847544 

.8225176822 

.8322216683 

2.3678227 

135.666296 

1.87224724 

-.8191344776 

.9285172734 

.3288559662 

2.32133516 

133.802755 

£ 

1.17681197 

-.8158857794 

.8186945753 

.8245325322 

2.27514727 

138.354333 

1.29154968 

-.8131886531 

.8170338888 

.82154277 

2.22964625 

127.749346 

1 

1.41747418 

-.8189494515 

.8155285649 

.8189941681 

2.18519888 

125.202259 

» 

1.55567417 

-9.89842782E-93 

.0141417407 

.9168114621 

2.14299359 

122.733252 

1.78735267 

-7.5479322  IE-83 

.8128854456 

.8149329382 

2.10863635 

129.35764 

»J5 

1.37381745 

-6.24547872E-83 

.8117487381 

.91338803 

2.86191571 

113.387544 

[•» 

2.85651234 

-5.28137654E-83 

.0186977233 

.311395411 

2.82339181 

115.931852 

a 

2.25781976 

-4.31843929E-93 

9. 7473470 4E-83 

.9184612495 

1.98794581 

113.394362 

2.4778764 

-3.58545172E-83 

8.88143782E-03 

9.57785998E-83 

1.95448981 

111.984057 

2.71858829 

-2.9767858  IE-93 

3.09243536E-03 

3.62254583E-83 

1.92327327 

118.195431 

3 

2.98364729 

-2.47131279E-43 

7.3735257E-03 

7.77664891E-03 

1.39419826 

108.529144 

3.27454922 

-2.85172729E-93 

6.71848199E-03 

7.82478344E-83 

1.36718448 

186.981955 

i, ,  < 

3.59381373 

-1.78338884E-83 

4.12163857E-03 

6.35428053E-83 

1.34213788 

185.549532 

k* 

3.94428413 

-1.41417483E-83 

5.57780 178E-03 

5.75428244E-83 

1.31918068 

184.226829 

c. « 

4.32376137 

-1.17487378E-83 

5.88228524E-83 

5.21413578E-83 

1.79782662 

183.887915 

C-4 

4.75881826 

-9.74736682E-84 

4.63878902E-03 

4.73224358E-03 

1.77825322 

101.386727 

a>! 

5.2148884 

-8.89243516E-84 

4.21948249E-03 

4.29430451E-83 

1.76828637 

108.357814 

5.72236778 

-4.71848184E-84 

3.8445624E-03 

3.90282442E-03 

1.74389221 

99.9125428 

*.v 

4.23829158 

-5.57788 164E-04 

3.58382287E-93 

3.5471513E-93 

1.72869999 

99.847197 

:;0 

6.89261226 

-4.6387889  IE-84 

3.19182324E-03 

3.22524072E-03 

1.71437394 

98.2558754 

*«*• 

7.56463344 

-3.84454231E-04 

2.98827846E-83 

2.93357183E-03 

1.78222342 

97.5305389 

*  *■<* 

3.38217587 

-3.19182317E-04 

2.44990773E-03 

2.44984132E-03 

1.49044924 

94.3632463 

m 

9.11142777 

-2.64998747E-04 

2.41449723E-03 

2.42899585E-03 

1.48819876 

94.2631754 

18.8888992 

-2.19999992E-04 

2. 19999995E-03 

2.21097261E-03 

1.67846593 

95.71863 

<• 

.  - 

Figure  19 
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AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  +-6Q  DECIBELS 
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Figure  20  Compensated  System  Frequency  Response 
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PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 


RAOIAN  FREQ. 

REAL  PART 

IMAGINARY  PART 

MAGNITUDE 

PHASE(RAD) 

PHASE ICEG) 

. 184741575 

-2.39220771 

-1.66091537 

2.91226675 

-2.53470893 

-145.2277:7 

.1149757 

-2.05136854 

-1.50709532 

2.54547622 

-2.50797801 

-143.494148 

.126185688 

-1.76821783 

-1.36634283 

2.23461116 

-2.48370669 

-142.335942 

.138488637 

-1.5329249 

-1.23744316 

1.97005693 

-2.46245877 

-141.333887 

.151991109 

-1.33732101 

-1,11928564 

1.74391164 

-2.44471727 

-140.372832 

.166810054 

-1.17461586 

-1.01885445 

1.54969324 

-2.43898665 

-139.235325 

.183073829 

-1.03916202 

-.911220404 

1.38289273 

-2.42149891 

-138.753176 

.200923381 

-.926259056 

-.819533617 

1.23676643 

-2.41725149 

-133.498358 

.220513875 

-.331990732 

-.735017078 

1.11016156 

-2.41800032 

-133.541263 

.242012828 

-.7530895 

-.656961126 

.999370661 

-2.42424292 

-138.900083 

.26560878 

-.686823804 

-.584718837 

.902018563 

-2.43432243 

-139.591057 

.291505308 

-.630904245 

-.517702353 

.316122474 

-2.45443313 

-148.62371 

.319926716 

-.533405553 

-.455388223 

.740089985 

-2.47882133 

-142.326051 

.351119176 

-.542701741 

-.397275798 

.672572108 

-2.50948973 

-143.794481 

.385352862 

-.507412466 

-.342966767 

.412448866 

-2.54721596 

-145.944774 

.422924291 

-.476358936 

-.292085857 

.558777222 

-2.59155843 

-148.434955 

.464158887 

-.44852822 

-.244322736 

.518755483 

-2.64231151 

-151.421999 

.509413806 

-.42304507 

-.199427017 

.46749463 

-2.78187883 

-154.760472 

.559081023 

-.399150742 

-.157212169 

.428995314 

-2.74638484 

-158.502233 

.613590733 

-.376188516 

-.117559888 

.394129591 

-2.33878547 

-162.645981 

.673415072 

-.353595756 

-.0804242261 

.34262655 

-2.91795858 

-147.134313 

.739072211 

-.330902379 

-.0458343711 

.334041632 

-3.88395521 

-172.114017 

.81113834 

-.387735349 

-.0138945796 

.388048867 

-3.89447228 

-177.414856 

.890215095 

-.293823366 

.0152206899 

.234234184 

3.08301791 

176.938457 

.977009969 

-.259835071 

.0412779193 

.242383326 

2.98356871 

178.945954 

1.07226724 

-.233343053 

.0640054395 

.241942139 

2.37387 962 

144.641232 

1.17681197 

-.20688473 

.8831222496 

.222958741 

2.75955236 

158.11876 

1.29154968 

-.179940328 

.0983768522 

.20507688 

2.44127174 

151.333777 

1.41747418 

-.152928819 

.189591363 

.188141559 

2.51979887 

144.373892 

1.55567617 

-.12637755 

.116795736 

.172021842 

2.39596143 

137.273533 

1.70735267 

-.108886423 

.119814167 

.15663175 

2.27064179 

138.898238 

1.37381745 

-.0770618877 

.119135249 

.141928356 

2.14475823 

122.885439 

2.85651234 

-.8554563234 

.115259292 

.127996638 

2.81924574 

115.694381 

2.25701976 

-.936506795 

.188620508 

.114591277 

1.89583025 

108.577274 

2.4770764 

-.0204899108 

.0999475468 

.102026215 

1.77308134 

101.58553 

2.71858829 

-7. 500239 15E-03 

.8899511678 

.0902633158 

1.65393513 

94.7664037 

2.98364729 

2.54491874E-03 

.0793108383 

.8793514584 

1.53871942 

38.1621603 

3.27454922 

9.37769306E-03 

.0686214773 

.9493287529 

1.4273337 

31.303874 

3.59381373 

.8148366908 

.9583582695 

.0682147409 

1.32183415 

75.7354599 

3.94420613 

.0178187351 

.0488619345 

.0520895757 

1.22110813 

69.9643471 

4.32876137 

.0192341152 

.0483426486 

.8446931816 

1.12590582 

44.5094745 

4.75081026 

.8194712689 

.0328972762 

.938227753 

1.83436815 

59.3795424 

5.2140084 

.8188733774 

.926533782 

.0325614184 

.952529207 

54.575923 

5.72236778 

.0177264567 

.9211973958 

.0276324635 

.374333383 

58.3954304 

6.28829158 

.8162567436 

.0167939828 

.8233734881 

.801451719 

45.9312765 

6.89261226 

.0146344848 

.9132103683 

.0197150191 

.734297995 

42.8721911 

7.56463344 

.8129813711 

.8103277974 

.9165885321 

.472043691 

38.5052889 

3.38217587 

.0113794273 

8.03199572E-03 

.0139285434 

.41463121 

35.2157369 

9.11162777 

9.87985853E-03 

6.21876925E-03 

.0114741036 

.541735109 

32.1879273 

10.0000002 

8.51098595E-03 

4.79675383E-03 

9.74963303E-03 

.5132213 

29.485425 

Figure  20 
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Figure  20  (Cont.) 


PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

XX  BODE  PLOT  (AMPLITUDE)  XX 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  COMMON  LOG  OF  AMPLITUDE 

MINIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  -  10  RADIANS/SEC 
AMPLITUDE  LIMITS  OF  BODE  PLOT  ARE  *-60  DECIBELS 


Figure  20 


ABSCISSA  ->  COMMON  LOG  OF  FREQUENCY 
ORDINATE  ->  PHASE  (DEGREES) 

TIC  MARKS  SHOW  MULTIPLES  OF  98  DEGREES 
MINIMUM  FREQUENCY  SHOWN  (X  ABSCISSA  -  .1  RADIANS/SEC 
MAXIMUM  FREQUENCY  SHOWN  ON  ABSCISSA  «  18  RADIANS/SEC 
MAXIMUM  PHASE  ON  ORDINATE  SCALE  -  180  DEGREES 
MINIMUM  PHASE  ON  ORDINATE  SCALE  -  -188  DEGREES 


Figure  20 
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ROOT  LOCUS 

PROBLEM  IDENTIFICATION  -  KATZ  EXAMPLE 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S  OR  W ' 


1 

? 


1 


OPEN-LOOP  ZEROS 
REAL  PART 

-.  1 
1 

DENOMINATOR  COEFFICIENTS 

0 

0 

4.6225 

4.3 

1 


IMAGINARY  PART 

0 

0 

IN  ASCENDING  POWERS  OF  S  OP.  W' 


OPEN-LOOP  POLES 

REAL  PART  IMAGINARY  PART 


-2.15  0 

-2.15  0 

0  0 

0  0 


MIN  GAIN 
0 


MAX  GAIN 

-30 


OPTION  TAKEN 

SIGMA  MAX  =  1 
OMEGA  MAX  =  1 


SIGMA  MIN  ■  -1 
OMEGA  MIN  *  -1 


Figure  21  Compensated  System  Root  Locus 
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ABSCISSA  ->  REAL  < SIGMA)  AXIS 
ORDINATE  ->  IMAGINARY  (OMEGA)  AXIS 
TIC  MARKS  SHOW  INTERVALS  OF  l 
THE  PLOT  FRAME  LIMITS  ARE: 
ABSCISSA,  -8  TO  2 
ORDINATE ,  -5  TO  3 


Figure  21 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  demonstrates  that  common  classical  control 
programs  may  be  adapted  to  run  on  an  inexpensive  microcomputer 
system.  With  the  exception  of  the  memory  management  problem 
encountered,  this  was  done  with  relative  ease.  This  allows 
more  people  access  to  these  programs  by  introducing  a  new 
group  of  computers  on  which  they  can  be  run.  Also  these 
programs  require  no  knowledge  of  any  computer  language  or 
input  card  formats  making  them  easier  to  use  than  the 
Fortran  versions. 

It  is  also  demonstrated  that  the  transfer  function  pro¬ 
grams  are  a  useful  tool  in  the  study  of  sampled  data  systems 
as  well  as  classical  systems. 

Only  a  sampling  of  five  programs  was  converted  in  this 
thesis.  Other  programs  can  and  should  be  similarly  converted 
to  run  on  microcomputer  systems.  There  are  programs  in 
existence  for  the  same  microcomputer  system  used  in  this 
thesis  that  generate  the  aircraft  stability  derivatives. 

An  effort  should  be  made  to  modify  these  programs  and  the 
programs  of  this  thesis  to  make  them  compatible  and 
complementary . 


APPENDIX  A 


DESCRIPTION  OF  MICROCOMPUTER  SYSTEM 

The  microcomputer  system  used  in  developing  this  thesis 
consisted  of  the  following  components: 

Apple  II  plus  computer  (48K) 

Disk  II  5  1/4"  floppy  disk  drive  and  controller  card 

USI  9"  green  screen  monitor 

NEC  PC-8023A-C  dot  matrix  printer 

Grappler  printer  interface  card 

Add  Ram  16K  expansion  card 

The  programming  language  used  was  Applesoft  basic.  All 
graphs  were  generated  using  the  High  Resolution  graphics 
commands.  In  the  High  Resolution  graphics  mode  a 
matrix  of  dots  280  dots  wide  and  192  dots  high  can  be  dis¬ 
played.  The  High  Resolution  page  one  is  the  only  page  used 
in  this  thesis  and  it  resides  in  memory  in  the  8,192-byte  area 
from  $2000  to  $3FFF.  High  Resolution  page  two  resides  in 
the  area  from  $4000  to  $6000. 

A  memory  management  problem  was  encountered  during  the 
programming  of  the  programs  requiring  the  use  of  the  graphics 
capabilities  of  the  microcomputer  system.  To  understand  the 
nature  of  the  problem  encountered  a  brief  description  of  the 
normal  use  of  memory  in  the  Apple  II  is  necessary.  An  Apple¬ 
soft  program  is  normally  loaded  at  memory  location  $800  and 
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loads  up.  LOMEM  is  set  to  the  end  of  the  program.  Simple 
variables  are  stored  from  LOMEM  upward  as  they  are  defined 
in  the  program.  This  gives  only  6K-bytes  of  RAM  before  there 
is  a  conflict  with  the  first  High-Resolution  page  and  only  14K 
before  there  is  a  conflict  with  the  second  High-Resolution 
page.  [Ref.  8]  Since  the  programs  and  variables  in  this 
thesis  exceed  14K,  simply  using  High-Resolution  page  two  in¬ 
stead  of  one  is  not  the  answer.  The  first  part  of  the  solu¬ 
tion  was  to  get  the  Applesoft  programs  to  load  above  the 
space  used  by  High-Resolution  page  one  and  in  effect  pro¬ 
tecting  that  space  from  interference.  This  fix  creates  a 
new  problem.  The  disk  operating  system  loads  at  the  top  of 
the  48K  of  memory  and  sets  HIMEM  to  $9CF8.  String  variables 
start  at  HIMEM  and  build  down.  This  allows  22K  of  space  for 
the  program  and  all  variables.  Some  programs  in  this  thesis 
require  more  space  than  this.  Even  with  a  16K  memory 
expansion  card  installed  the  disk  operating  system  will 
ignore  it  so  the  added  memory  is  useless.  This  final  problem 
was  solved  by  using  a  utility  program  that  relocated  the  disk 
operating  system  into  the  higher  memory  provided  by  the 
memory  expansion  card  and  resets  HIMEM  to  SBFOO.  This  gives 
32. 5K  of  useable  program  space  which  is  sufficient  for  all 
programs  in  this  thesis.  [Ref.  9] 
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